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FOREWORD 


The research presented in this report was performed for the Astrionics 
Laboratory of the George C. Marshall Flight Center, Huntsville, Alabama. 
This report is the final report on the interplanetary navigation and 
guidance study task under NASA Contract NAS 8-20358. The results pre- 
sented here extend the scope of the navigation and guidance study completed 
under Contract NAS 8-11198. 
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ABSTRACT 

The primary objective of the Navigation and Guidance research task of 
Contract NAS 8-20358 is to critically assess the influence of the constraints 
imposed on the study performed under Contract NAS 8-11198 on the results 
obtained. The study has analyzed the influence of bias error sources on the 
navigation and guidance requirements for a round-trip Mars mission. The 
equation of motion bias errors considered are: (1) Astronomical Unit Con- 

version, (2) Solar-Radiation Pressure, (3) Mars ' Planetary Mass , and 
(4) Earth's Planetary Mass. The measurement bias errors considered are: 

(1) Sextant Angle Measurement Bias, and (2) Onboard Clock Bias. The 
navigation system used in the analysis consists of a sextant with a 10 arc 
second accuracy. The measurement data are processed with a Kalman filter. 

The theoretical analysis of the effect of bias error sources on the estima- 
tion process is presented. Recent theoretical developments are presented 
that concern three areas of interest in this study. An analysis is pre- 
sented that shows the error introduced in an estimation process due to 
neglecting both equation of motion and measurement bias error sources in the 
modeling of the physical process. A second analysis shows the capability 
of separating the effects due to bias error sources from those due to ran- 
dom errors. Finally, techniques are presented for efficient parametric 
analysis by means of matrix manipulations. 

Data results are presented that show the effect of bias errors on the navi- 
gation and guidance requirements. The results are presented in terms of 
end point constraint estimate accuracies and velocity requirements, res- 
pectively. The only equation of motion bias error that significantly limits 
the navigation system performance is the uncertainty in the mass of Mars. 

An uncertainty of 150 km 3 /sec 2 , produces an uncertainty in the entry altitude 
at Mars of 4 to 5 km. The two measurement bias errors contribute significant 
estimation errors when they are neglected. Expanding the state being esti- 
mated to include the measurement bias error sources, allows the navigation 
system to be calibrated with the resulting improved performance. 
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The midcourse guidance velocity requirements are significantly altered by 
the equation of motion errors. The trajectory deviations due to the bias 
errors cannot be accurately estimated until the final portion of the tra- 
jectory at which time they must then be corrected by the guidance system. 
An example cited for an Earth-Mars trajectory shows an increase in the 
velocity requirements from 23 meters/second with no equation of motion 
errors to 41 meters/second with uncertainties in the A.U. conversion and 
Mars' planetary mass . 
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SECTION I 
INTRODUCTION 

1.1 GENERAL OBJECTIVES AND SCOPE 

The primary objective of the Navigation and Guidance research task of Con- 
tract NAS -8-20358 is to critically assess the influence of the constraints 
imposed on the study performed under Contract NAS-8-11198 on the results 
obtained^ . Research under Contract NAS-8-11198 established the basic 
requirements for an Advanced Spaceborne Detection, Tracking and Navigation 
System capable of performing future interplanetary missions. The constraints 
on the original study that are analyzed in this research report are the 
following: 

(1) The assumption that the only significant error sources are 
random in nature. 

(2) The original choice of 1975 round-trip Mars trajectory. 

In order to analyze the influence of these constraints on the previous 
results, several tasks had to be completed: 

(1) Research and analysis of the literature available on the 
techniques used to estimate the various physical constants 
that are required to specify an interplanetary trajectory. 

The uncertainty associated with each of these estimates con- 
stitutes a bias error in the equations of motion of a spacecraft. 

(2) Determine the types of measurement bias errors to be considered. 


( 1 ), 


Superscripts refer to references listed in Section 8. 
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(3) Derive suitable mathematical models to analyze the influence 
of the equation of motion and measurement bias errors on the 
trajectory estimation process. 

(4) Develop methods for data presentation which indicate the 
importance of the bias error sources. 

(5) Determine the areas which require future research. 

The scope of the study includes an evaluation of the onboard navigation 
system performance under the influence of the following bias errors. 

(1) Equation of Motion Errors 

a. Uncertainty in the Astronomical Unit conversion 

b. Uncertainty in the Earth's mass 

c. Uncertainty in Mars' mass 

d. Uncertainty in the solar radiation pressure 

(2) Measurement Biases 

a. Sextant bias in onboard measurement 

b. Time bias in onboard clock 

The results that are obtained show the effects of these bias errors on the 
navigation system performance. In order to make the problem amenable to 
study, certain restrictions on the scope had to be made. 

(1) The covariance matrix of injection errors at Earth is not 

studied as a parameter. This matrix, which is a function of 
the time in park orbit at Earth, is intended to be represen- 
tative of the capabilities of future launch vehicle guidance 
systems. The primary influence of this matrix is on the 
magnitude of the midcourse velocity requirements at the first 
guidance correction. 
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(2) The study emphasizes the following phases of the mission: 

a. Midcourse from Earth to Mars 

b. Midcourse from Mars to Earth. 

(3) The onboard navigation instrument is a sextant with a random 
measurement error of 10 arc seconds. The error magnitude is 
not studied parametrically. This random error is a parameter 
in the original study. The measurement schedules used are 
also derived from the original study. 

(4) The guidance system analysis is restricted to evaluation of 
trajectory deviation corrections rather than a statistical 
analysis. This is due to both theoretical and simulation 
limitations in the analysis of the effects of bias errors 
on the guidance system performance. 

1.2 STUDY FORMAT 

The effect of bias errors on the navigation system performance may vary 
considerably depending on the mission itself. This study is designed to 
identxfv the error sources and show their effect on the system performance 
for a variety of trajectories. 

The study was performed in four basic steps. The results of each step are 
presented in the order that they were performed. 

(1) Analysis of applicable filter theory 

(2) Analysis of bias error sources. 

(3) Computer simulation. 

(4) Generate results. 

Section 2 contains a description of the theoretical aspects of Kalman 
filtering with bias errors present in the system. There are several new 
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and interesting developments presented that are incorporated in the study. 
These developments concern the following areas of interest. 

(1) The effects of neglecting bias errors. 

(2) Separable properties of bias error effects from the effects 
of random errors. 

(3) Computer end of run observations that greatly s imp lify 
parametric type analyses. 

Section 3 presents a description of the bias errors that are considered 
in the study. A particular emphasis is placed on the equations of motion 
parameters. The uncertainty in the astronomical unit conversion is dis- 
cussed in detail and a model described for analyzing its effect on an 
interplanetary trajectory. 

The digital computer simulations that are utilized in the study are des- 
cribed in Section 4. 

The results of the parametric analyses for the outbound and return trajec- 
tories are presented in Section 5. 

The tables and figures that are discussed in the report are presented at 
the end of the report. 
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SECTION 2 

THEORETICAL ANALYSIS OF OPTIMAL ESTIMATION PROCESS 


The function of the navigation system, as defined in this report, is to 
obtain an estimate of vehicle state based on either direct observations of 
the vehicle (Earth-based tracking) or observations of celestial bodies 
whose positions are known (on-board tracking). The purpose of this section 
is to describe the error analysis techniques that have been used to study 
the effects measurement and equation of motion parameter uncertainties on 
the vehicle state estimation process. 


The following presentation of the orbit estimation process is based on the 

( 2 ) 

filtering theory of Kalman and the orbit determination application by 
q » . j . (3)4)5) 

bcnmiut . The theory for including measurement and equation of 

motion bias type errors into the optimal orbit estimation process is pre- 
sented in References 1 and 6. These references show three different 
techniques by which the bias errors in the estimation process can be handled. 

technique produces an optimal estimate but under a different assump- 
tion concerning the modeling of the dynamic and/or measurement bias error 
sources. The three techniques are the following. 


O) Neglect - In this case, the state being estimated is simply 
the three components of vehicle position and velocity. The 
dynamic and measurement models are assumed to be known per- 
fectly and there are no bias errors in the physical process 
being analyzed. The effect of uncertainties in system param- 
eters are neglected . 

(2) Include - This technique allows for an expansion of the state 
vector to include specified dynamic and measurement parameters 
that produce bias error sources. The optimal estimation 
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process fits the measurement data to the model to obtain a best 
estimate of the vehicle state and the system parameters 
simultaneously. 

(3) Consider - With this type of analysis, the measurement and 
dynamic parameters are not included in the state vector being 
estimated but the effects of the uncertainties in these param- 
eters are considered in the estimation process. An optimal 
estimate of the state is obtained under the constraint that 
specified unsolved for system parameters have constant uncer- 
tainties (standard deviations) associated with their fixed 
estimates. 

Each of these techniques has both favorable and unfavorable aspects. If 
the bias errors are neglected, the state vector is reduced in dimension 
and correspondingly the number of computations required in the estimation 
process is reduced. Alternately, the removal of bias errors from a problem 
in which they can significantly influence the results being obtained causes 
the state estimates that are obtained to have optimistic uncertainties 
attached to them. The implication here is that an analysis must be made 
to establish the relative importance of measurement and dynamic bias error 
sources to a particular orbit estimation process. On the basis of this 
analysis, the method for treating each of the error sources can be esta- 
blished. The allowable total size of the state vector would be limited by 
the computer capability either ground or onboard in terms of its size and 
speed. Within these limitations, the filtering techniques used for the 
various bias errors sources in the order of decreasing importance would be 
to include, consider, and neglect them. 

The second part of this section describes some recent developments on 
techniques for efficient parametric analysis of navigation system per- 
formance. The theory is based primarily on the fine work of Gunckel^. 

The techniques apply to two areas: the effect of bias errors that are 
not included in the optimal filter, and the effect of varying the initial 
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state variable (bias error) uncertainty covariance matrix. The important 
| feature of these techniques is that they are applied to the data obtained 

at the conclusion of a single computer run of an orbit determination 
I simulation. The parametric analysis is accomplished by means of matrix 

manipulations of the end point data and thereby eliminates the need for 
• making a large number of computer analysis runs when scanning a range of 

parametric values of the initial error variances. 

■ 2.1 SYSTEM DESCRIPTION 

| The Kalman filter optimal estimation process is based on minimizing a 

loss function for a linear dynamical system. In order to get a linear 
I system of equations for the orbit determination process, a linearization 

is performed about the "best estimate" of the state of the nonlinear system. 

The equations of motion of the spacecraft may in general be written as a 

■ set of nonlinear differential equations of the form 

| X = F(X, U, T) (2-1) 

I The n dimensional vector, X, defines the vehicle state, u is an f dimen- 

sional vector of control and forcing functions. Uncertainties in the 

■ functions, U, will be referred to as "equation of motion bias errors". 

The observations or measurements are in general related to the state X by 

Y = G(X, V, t) + q(t) (2-2) 

a The m dimensional vector, V, defines the measurement parameters and 

■ q is the random measurement noise. Uncertainties in the vector, V, will 
be referred to as "measurement bias errors." Linearization of equations 

S 2-1 and 2-2 about a nominal trajectory yields the following. 
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The solution of (2-3), for constant 
may be expressed as 

x(t 2 ) = cp(t 2 , tj) x ( 


A(t) x + B(t) u (2-3) 

s H(t) x + D(t) v + q(t) (2-4) 

controls over the interval t^t^t 2 , 

+ W t l ) u (t l ) (2-5) 


where cp is the nxn state transition matrix that relates a deviation in the 
state at t^ to a deviation at t 2 » and cf> u is an nxA matrix that relates the 
deviation in the state at t 2 to a unit variation in the control at t.. 

The elements of 9 and are referred to as sensitivity coefficients. 

In general 9 and cp^ are computed by numerical integration of the varia- 
tional equations with appropriate initial conditions. 

The solution of the linearized differential equations of spacecraft motion 
(equation 2-5) and the state measurements or observations (equation 2-4) re- 
present the physical process shown in Figure 2-1. The figure is drawn 
with fat lines to represent the fact it is a matrix block diagram. The 
three procedures that are used to model and estimate the bias errors shown 
in Figure 2-1 are discussed in the next sections. 

2.2 KALMAN FILTER THEORY 

Derivations of the optimal Kalman filter for different dynamic and measure- 
ment models are presented in many references. The theory presented 
in the following paragraphs will draw heavily from techniques used in 
References 6, 7, and 8. 
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2.2.1 Bias Errors Neglected 

The navigation analysis performed during the initial study ^ neglected 
the effects of the bias errors, u and v, shown in Figure 2-1. The optimal 

/ Q\ 

filter for the case of neglecting these errors will be derived below v . 

The derivation will be followed by an analysis^ Q f the error in the 
estimate due to the fact the bias error sources are neglected in the filter. 

The recursive linear estimation problem may be stated as follows. Given 

the physical process shown in Figure 2-1 determine an estimate, x (t) , of 

n 

the state at t that is a linear combination of an estimate at t-1 and 
the measurement data y. The estimate must be "best" in the sense that 
the expected value of the sum of the squares of the error in the estimate 
is a minimum (minimum variance estimate). That is, x^(t) is to be 
chosen so that 


E 




minimum 


( 2 - 6 ) 


The dynamics of the system are described by the homogeneous linear difference 
equation 


x(t) - co(t, t-1) x(t-l) (2-7) 

The measurements, y, are linearly related to the state and corrupted by 
additive noise. 


y(t) = H(t) x(t) + q(t) 


(2-8) 


The Kalman filter equation shall be derived for the model described by 
equations (2-7) and (2-8). 


* e[ ] means expected value of bracketed quantity. 
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The form of a linear estimation equation can be hypothesized from the 
physical characteristics of the system. The state evolves according to 
(2-7) so, given an estimate x n at t-1, it is reasonable to predict the 
estimate at t as 

st(t) = cp(t, t-1) * n (t-l) (2-9) 

when no information is available. A measurement at t can be used to 
modify the estimate. Based on x(t) and (2-8), one would expect the measure- 
ment value at t to be H(t) x(t). An error in the estimate is reflected 
by an error in this expected measurement value. 

e(t) = y(t) - H(t) cp(t, t-1) 3^(t-l) (2-10) 


The estimate is to be a linear function of the new measurements. Define 

an unknown "gain" matrix, K(t), such that the new estimate 5c at t after 

n 

inclusion of the observation is given by 


s n (t) - 


— 

3 


x(t) + K(t) y(t) - H(t) S(t) 


■■ 

) 


( 2 - 11 ) 


The linear filter described by equations (2-9) and (2-11) is shown in 
block form in Figure 2-2. 


The matrix K(t) shall be determined so that E 


[jx n (t) - x(t)j T ^ n (t) - x(t)j 


is 


minimized. It shall be referred to as the weighting or gain matrix. 


Let 


x (t) = x (t) - x(t) 

n n 


(2-12) 
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then 


E [ V t> - x(t)) T (^(t) - x<t>) • E x n (t) T x n (t) 


(2-13) 


This can be rewritten as 


E [*n^ T “ trace E^Ct) x^(t)^j 


(2-14) 


where the trace is defined as the sum of the diagonal elements of a matrix. 


Define the matrix P (t) as 

n 


P n (t) - E 


[yt) ?co T J 


(2-15) 


Now, form x n (t) using equation (2-11) 


^ n (t) = |tp(t,t-l) x(t-l) + K(t)(y(t) - H(t) cp(t,t-l) x(t-]))j - ro(t,t-l) x(t-l) 
= <p(t,t-l) x(t-l) - K(t) H(t) cp(t, t-1) x(t-l) + K(t) ^H(t) x(t) + q(tj) 
= (l - K(t) H(t)) cp(t, t-1) x(t-l) + K(t) q(t) (2-16) 


From equation (2-16) the P (t) matrix can now be formed 

n 


?„(t) - e| 


^ j<I-i 


(I - KH) o(t, t-1) x n (t-l) + Kq| |(I-KH) cp(t,t-l) X n (t-1) + Kq 


I 


- (I - KH) fp(t 


+ K E 


,t-l) E^ n (t-1) x n (-l) T ] t " 1 > " hT rT > 

|q * n (t-D T J w T (t, t-1) (I - H T K T ) 

+ (I - KH) cp(t, t-1) E^ n (t-1) q*J K T + K E |q q T JlC T 


(2-17) 
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By definition 


q(t) q (t) - Q(t) 


x (t-1) x 
n 


q(t) q (t 


ejq(t) ^(t-i) 


i 

-] 




0 = E 


t(t-D q(t ) T ] 


Therefore P (t) can be rewritten 
n 


P n (t) = (I - KH) P(t) (I - KH) T + KQK T 


where 


(2-18) 


P(t) " cp(t, t-1) P n (t-1) cp A ( t , t-1) 


(2-18a) 


Expanding (2-18) gives 

P R (t) = P(t) - KHP(t) - P(t)H T K T + K(HP(t) H T + Q) R T (2-19) 

The matrix P(t) doe* -rot depend upon K(t), so it is unaffected by the 
selection of K(t). The matrix (HP(t) H + Q) is symmetric and nonnegative- 
definite, so it can be written as the product of a matrix S(t) and its 
transpose (i.e., a matrix square root). 

S(t)S(t) T “ HP(t)H T + Q (2-20) 

The last three terms of (2-19) have the form of a quadratic matrix poly- 
nominal in terms of the unknown K(t). Introduce (2-20) into (2-19) and 
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hypothesize the existence of a matrix R(t) such that 

P n (t) « P(t) + (KS - R) (KS - R) T - RR fc (2-21) 

This procedure is the matrix equivalent of completing the square of a 
quadratic polynominal. Assuming that S(t)S(t) is positive definite, it 
follows directly that 


R(t) - P(t) H T |s(t) _1 j T (2-22) 

Only the product term in (2-21) involves the gain matrix K(t). The 
product of a matrix and its transpose is nonnegative-definite, so the trace 
of P n (t) is minimized by choosing 


K(t)S(t) 


P(t) H(t) T 



Thus the optimal gain matrix is 


(2-23) 


K(t) 


P(t) H x (t) 


H(t) P(t) 


T 

H A (t) 


+ Q(t) 


3 


-1 


(2-24) 


Substituting (2-24) and (2-22) into (2-21) shows that the convariance 
matrix P n (t) for the optimal gain is 


V c) 


- poo - p(t> 


T 

H (t) 


H(t) P(t) 


T 

H (t) 


i 


-1 


+ Q H(t) P(t) (2-25) 


Equations (2-9), (2-11), (2-18a), (2-18), and (2-24) constitute the Kalman 
filter for the model described by Equations (2-7) and (2-8). These equa- 
tions are summarized below. 
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The time propagation of the state estimate is given by equation (2-9). 


x(t) ■ cp(t.t-l) x n (t-l) 

The updating of the state estimate for a measurement is given by equation 

( 2 - 11 ). 


x n (t) 


= x(t) + K(t)[y(t) - H(t) x(t) 


L. 


— 

) 

— 


The time propagation of the covariance matrix is given by equation (2-18a). 


P(t) - w(t, t-1) P n (t-1) cp(t, t-l) T 

The updating of the covariance matrix for a measurement is given by 
equation (2-18). 


P n (t) = (I - KH) P(t) (I - KH) T + KQK T 


The optimal filter gain is given by equation (2-24) 


1 


K(t) = P(t) H(t) x H(t) P(t) H(t) r + Q(t) 


-1 


2 .2. 1. 1 Effect of Neglecting Bias Errors . The bias errors, u and v, 
shown in Figure 2-1 were neglected in the derivation of the optimal filter 
in the previous section. This section will present an analysis, due to 
Gunkel^, of the error in estimate caused by neglecting these bias errors. 


Two different types of bias error sources are shown in Figure 2-1, state 
vector perturbation biases, u, and measurement biases, v. State vector 
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biases might include uncertainties in planetary gravitational constants 
and solar pressure acceleration. Measurement biases could include any of 
a number of instrument imperfections. 

The analysis of the effect of these error sources is acconplished by an 
application of linear superposition. The system and optimal estimation 
process models are linear dynamical systems. Therefore, the effect of any 
additional inputs can be found by ignoring the effects of any of the original 
inputs (the random noise q(t)). 

Equation of Motion Biases . The effect of the equation of motion biases 
is shown in the following equations 


x u (t)= w(t, t-1) x^(t-l) + r (t, t-l)u (2-26) 

*nu (t) = I" 1 " K(t)H(t)j cp(t, t-1) x.^t-1) + K(t)H(t) x^t) (2-27) 


where the subscript U> is used to identify the value of the state vector 
resulting from the input m>. 

Equations (2-26) and (2-27) can be obtained either by examining Figures 2-1 
and 2-2 or by using the model equations given in the previous section. 

The estimation error resulting from u is found by taking the difference 
between (2-27) and (2-26) to give 


*nu (t > * V° ‘ X - (t) = 


I - K(t)H(t) 


n(t,t-l) x (t-1) - co (t, t-l)u 

“u u< 


(2-28) 
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One of the primary objectives of this study, is to determine the sensitivity 
of the estimation error to changes in the state vector biases. These sensi- 
tivities are found by partial differentiation of (2-28). Defining the sen- 
* 

sitivity matrix, C (t), as 


V fc > 





and then differentiating (2-28) gives 


yt) = j^I - K(t)H(tj |o(t, t-1) C u (t-1) - cpjt, t-lj 


(2-29) 


The bias errors being analyzed are generally the result of an error in 
estimate of one or more physical constants in the equations of motion. 
The expected value of the error in estimate of the constants is zero. 


E(u) » 0 


(2-30) 


Due to the uncertainty in the estimates being used for the constants, u, 
there is a covariance matrix associated with the estimate. 


E(uu T ) = M (2-31) 

The uncertainties in the equation of motion constants (2-31) can be trans- 
formed into the covariance matrix of error in estimate of the state by 
means of (2-29). 


*The matrix has been defined as C(t) because an interesting relationship 
exists between this sensitivity matrix and the correlation matrix between 
the vehicle state and the system biases. This relationship will be shown 
in Section 2.2.3. 
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P u (t) " C u (t) M C u (t)T 


(2-32) 


where: P u (t) the covariance matrix of the error in estimate of the 

state due to the uncertainties in the equations of motion that have been 
neglected. 

Measurement Bias Errors . The effects of measurement biases can be deter- 
mined in a similar manner. The dynamical equations are 


y v (t) = D(t) v 


x (t) = 

nv 


•3 


I - K(t)H(t)| CJ(t, t-1) x nv (t-l) +K(t)D(t)v 


(2-33) 

(2-34) 


The sensitivity C v (t) is 


asr (t) 

c (t) = — — 

v x *v 


I - K(t)H(t) 


’(t, t-1) C v (t-1) + K(t)D(t) 


(2-35) 


As in the case of equation of motion errors, the expected value of v is 
zero but there is an associated error in this estimate. 

E(v) = 0 E(w T ) " W (2-36) 

The uncertainty in the estimate of the measurement biases, W, can be 
reflected in an error in estimate of the state with the use of (2-35). 


P (t) - C (t) W C (t) T 

V V V 


(2-37) 
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where P y (t) is the covariance matrix of the error in estimate of the state 
due to neglecting the measurement bias uncertainties. 

The total error in estimate of the state must realistically include the 
effects shown by equations (2-32) and (2-37) in addition to that shown by 
(2-25) which is the result of initial state uncertainty and random noise 
in the observations. 


P (t) = P(t) + p (t) + P (t) 

i u v 


(2-38) 


where P fc (t) is the total state uncertainty including the effects of neglected 
bias errors. 

Portions of the data presented in later sections were obtained by the tech- 
niques just presented. These data show the effects of the bias errors that 
had been neglected in the previous study^^. 

2.2.2 Bias Errors Included (Parameter Estimation ) 

The bias error parameters, u and v, can be included as additional ele- 
ments of the state vector and estimated along with the vehicle state. The 
physical process for the inclusion of the bias errors in the state is shown 
in Figure 2-3. 

The dynamics equations for the system may be written in terms of partitioned 
matrices as 



/®(t, t-1) co (t, 
i u 

0 I 


V 


o 


0 
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or by definition 


z(t) = cp 2 (t,t-l)z(t-l) 


(2-39a) 


and 


y z (t) - ( H(t) 0 D(t)J / x(t)\ + q(t) (2-40) 


and by definition 


u 


y z ( t) - H z (t) Z(t) + q(t) 


(2-40a) 


These equations have the same form as (2-7) and (2-8) and yield the same 
solution for the optimal filter as shown in Figure 2-2. 

The filter equations and optimal gain are presented below. 


N> 

/-N 

ft 

II 

cp z (t, t-1) z„(t-l) 

(2-41) 

P (t) = 

Z 

t-1) P zn (t-1) q£(t, t-1) 

(2-42) 

S n (t) = 

z(t) + K z (y z - y z ) 

(2-43) 

P zn (t > " 

P z (t) - K z H z P z (t) 

(2-44) 
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where the optimal gain, K (t) is 

z 

K (t) = P H T fH p H T + Q) 

2 2 Z K 2 2 Z 


(2-45) 


Equations (2-41) and (2-42) are for updating the state and covariance 

matrix between the observations, (2-43) and (2-44) are for the improvement 

in estimate and the covariance matrix of the error in estimate as a result 

of the observation y (t). 

z 

With this expanded state definition one may include as many unknowns (in 
principle) as he desires. The expanded state being estimated increases 
the computer storage requirements and the time required to perform the 
computation. 

2.2.3 Bias Error Considered 

The bias error treatment described in this section is due to Schmidt . 

An optimum estimation process is desired that includes the effects of bias 
error parameters on the state estimation process but constrains these 
parameter estimates and their uncertainties to remain constant. 

From the previous derivation, it is apparent that no theoretical difficul- 
ties are introduced by parameter estimation. There may be practical diffi- 
culties since one may not be able to obtain a solution that converges when 
a large number of unknowns are introduced. This can occur when the unknowns 
are not linearly independent for the number of significant figures retained 
in the numerical calculations. Also the computer size can be excessive when 
a large number of unknown parameters are added. 

For these reasons one would like to include the effects of unknown parameters 
in the sense that they deteriorate the estimate of state, without actually 
carrying through all the calculations for estimating them. The manner in 
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which this is done will be described subsequently for: (1) equation of 

motion bias errors and (2) measurement biases. 

The analysis will be performed independently for the two bias error sources 
with the random noise, q(t), neglected in both cases. This type of 
analysis invokes the principle of superposition in a linear system. The 
following analysis will be directed at three areas of interest: (1) obtain- 

ing an expression for the optimum gain, K(t), when the bias error effects are 
considered, (2) showing the separability of the total error covariance matrix 
into a part due to random errors and parts due to each of the bias error 
sources, and (3) show the relationship that exists between the correlation 
matrix of the state and bias errors and the partial sensitivity matrix 
obtained in Section 2.2. 1.1. 

The physical process being analyzed is shown in Figure 2-1 and the linear 
filter is shown in Figure 2-2. The filter gain, K(t), is to be selected 
so that it is optimum in the sense of a minimum variance filter with the 
effects of the bias errors, u and v, considered. The bias errors are 
assumed to be defined by the following statistical quantities. 


E(u) = u = 0 

E(v) s v : 0 


E(uu^) - M (2-46) 

E(w T ) ~ W (2-47) 


Equation of Motion Errors . The following equations describe the dynamics 
of the system and the estimation process with the equation of motion bias 
errors. 


x(t) 

*n (t > 


cp(t, t-1) x(t-l) + to (t, t-l)u 


u 




KH 


cp( t , t-1) x (t-1) + KH x(t) 
n 


■] 


(2-48) 

(2-49) 
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The optimum filter gain is to be selected so that E (x(t) 
V t))j is minimized. 




Let 


then 


* n (t) = " ft n (t) 


/ . r- 

(* (t) - x n (t)) T (x (t) - x n (t)J s E x n (t) T x n (t) 

L. — 


or 


: j5n <t)Ti » <C ] 


= Trace E 


* n (t) \, T(t >] 


” Trace 


M 


now form (t), 
n 1 


x n ( t) = 


I - KH 


I 0(t , t-1) x n (t-l) + r ^ u (t, t-l)uj 


From this the covariance matrix of the error in estimate, P (t) 

n 

formed as defined by (2-52). 


P n (t> - E, 


If ] 

it ' “J 


p(t, t-i) x n (t-i) + cr u (t, t-i)u 


|jp< t , t-l) x n (t-l) + v u (t, t-1) uj T - KhJ 

P n (t) - [l-KH]{ <Pe[ x n (t-l)x n T (t-l)j9 T + 9E[x n (t-l)u T ]q 
+ 9 u e[u x n T (t-l)] cp T + 9 u e[uu T ] q£ } [i - KhJ 
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(2-50) 


(2-51) 


(2-52) 


(2-53) 


can be 


(2-54) 

> T 

u 

(2-55) 
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At this point it is very useful to define several of the statistical 
correlation terms that are in (2-55). 

Let 


c u (t) ■ E |*(t) „ T ] c un (t) e[? o(£) jJ (2 . 56) 

using the dynamic equation (2-48) 

c u (t) = e| ^(t, t-1) x n (t-l) + ’? u (t, t-1) uJu T J (2-57) 
or using (2-56) and (2-46) 

C u (t) = ^(t, t-1) C un (t-1) + - u (t, t-1) M (2-58) 


The correlation factor, C (t), changes at the time of an observation. 

u 


C 

un 


(t) 



from (2-53) and (2-58) 


(2-59) 


C 

un 


(t) 



(2-60) 


The correlation matrix (2-60) is of the same form as the sensitivity 

matrix shown in equation (2-29). If equation (2-60) is post multiplied 

by M the two equations are the same. Interpretation of the correlation 

matrix as a sensitivity matrix suggests the following separated form for 

the covariance matrix P (t). 

n 
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let 


p n ( t -i) - p;<t-D + c un (t-i) m ’ 1 c^(t-D 


(2-61) 


where 


P (t-1) is the total state covariance matrix, 
n 

P^(t-l) is the uncertainty in the state due to the initial covariance 
matrix, P(o), and random measurement noise. 

The second term is the uncertainty in the state due to the uncertainty in 
the equation of motion parameters. 

Returning to equation (2-55) and making the substitutions for the expected 
values yields the following 


P n ( t ) - 


I - KH 


I'O P (t-1) cp T + c P C (t-1) cp T + cp c T (t-l)cp T +9 M cp T f [*I-KH 
In un T u u un u u J L J 


(2-62) 


Using the assumed form of P^Ct-l) in equation (2-61) for P (t-1) in 
(2-62) yields 


p„(t) 


* jl - KH 


;(t-l)cc T + ro C un (t-l)M _1 C un T (t-l)cp T 


+ cp C (t-1) r T + ro C T (t-1) cp T + w M *p T 
T un u u un ' ' T u u 


I - KH 


(2-63) 


Evaluation of C (t) M -1 C (t) T using (2-58) yields 


.-1 


•1 „T 


C (t) M C (t) = cp C (t-1) M C (t-1) co 1 + c (t-1) n + rc c 1 ( t _i) 
u u un' un un' / u u un' ' 

T 


ro 


+ cp M rp 

u u 


(2-64) 
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Substituting (2-64) into (2-63) yields 


V'> 


r 


I - KH 

P'(t) + C (t) m” 1 c T (t)\ 

r- • 

I - KH 

- 

L u u ] 



(2-65) 


where 


P'(t) = - P'(t-l) c? T 

n 


using (2-60) and (2-18) 


P n (0 * yt) + C un (t) M- 1 C^t) (2-66) 

Equation (2-66) which is the recursion of (2-61) shows the feasibility 
of expressing the total covariance matrix in a separated form. One part 
is due to the initial covariance matrix, P(0), and the random measure- 
ment noise and the other is due to the effect of the bias errors. 

The filter gain, K(t), is to be selected to minimize the trace of 
equation (2-65). The derivation of the solution for K(t) is the same 
as that presented in Section 2.2.1. 

Rewriting (2-65) 


P n (t) = (I - KH) A(I - KH) T (2-67) 


where 


A ^ P'(t) + C u (t) M' 1 cj(t) 
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or 


P fl (t) = A - KHA - AH T K T + KHAH T K T 


(2-68) 


This equation is exactly of the same form as (2-19) with the following 
equivalent variables. 


A ~ P(t) K ~ K H ~ H HAH 1 ~ SS T 


The solution of (2-68) for the optimal gain is given by (2-23). Making 
the appropriate substitutions into (2-23) yields 


K(t) S(t) 


A(t) H T (t) 



(2-69) 


or the optimal gain when considering the bias error parameters is 


K(t) = AH T (HAH T ) _1 


(2-70) 


where 


A = P'(t) + C u (t) M _1 C^(t) 

Coirparing the filter gain obtained while neglecting the bias errors 

(2-24) with (2-70), indicates that the only change in the gain matrix is 

due to a change in the form of the covariance matrix, P(t). In the gain 

computation for the inclusion of the bias error effects, the covariance 

matrix used is increased by the effect of the parameter uncertainties given by 
-IT 

the term C (t) M C (t). 

u u 

2-22 


PHILCO 


philco ford corporation 


Space Sl Re-entry 
Systems Division 



SRS-TR148 


Measurement Bias Errors . The method of considering measurement bias 
errors in the estimation process can be analyzed in a manner similar to 
that just described for the equation of motion biases. Linear superposi- 
tion will be invoked to neglect all inputs other than the measurement 
biases in this analysis. The physical process and filter being used are 
shown in Figures 2-1 and 2-2 respectively. 


The following equations describe the dynamics of the system and the esti- 
mation process with the measurement biases. 



x(t) 

■ t-1) x(t-l) 

(2-71) 


y(t) 

= D(t)v + H x(t) 

(2-72) 

p x> 

rt 

It 

"TT " 1 

1 

S 

rt 

rt 

t 

i—* 

x (t-1) + K |D(t) v + H x( t ) I 

(2-73) 

V- -J 


L J 



The error in estimate x n (t) is the following 


x n (t) = x(t) - X R (t) 


I - KhJ o(t. 


t’ 1 ) * n (t-l) - K D(t) v (2-74) 


The covariance matrix is 


P n (t > ‘ E 


x n (t) V fc > 


= E 


(I - KH) rc(t, t-1) x n (t-l) - K D(t) 


(I - KH) o(t, t-1) x* (t-1) - KD(t) v 

a 


■} 

(2-75) 
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yt) 


■ (I - KH) co(t, t-1) E 
- (I - KH) «(t, t-1) E 


x (t-1) x 

n 


x n (t-l) v T I d T K X 




cp T (t, t-1) (I - KH) T 


K D E 


V ,jT (t » 


t-1) (I - KH) + K 


D E 1 (v v T )J 


T T 
D K 


Let 


C (t) 
v 


:|x(t) v T 


c w (t > “ 


then from equation (2-71) 


C (t) 

v 


= E 


'P ( t , t-1) x (t-1) v J 
n 


o(t, t-1) C vn (t-1) 


and from equation (2-74) 


C (t) = E ( 

vn ‘ 




(I - KH) cp(t, t-1) x n (t-l) v 1 - K(t) D v 


applying (2-77) and (2-46) yield 


C vn (t) = (I " t_1 ) cyt-1) " K(t) D(t) W 


or 


C vn (t) = (I - KH) C v (t) - K(t) D(t) W 
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The covariance matrix will be defined to be made of two parts, one due 
to the random errors and one due to the biases. 

Define 

P n (t_1) = P n (t " 1) + c V n (t_1) W_1 (2-81) 


and 


P(t) = cp( t , t-1) P n (t-1) fn(t, t-1) 


(2-82) 


or using (2-77a) 


p(t) = p'(t) + c v (t) W ' 1 c*(t) 


(2-83) 


where 


p’(t) = 


'0(t, t-1) P^(t-l) t-1) 


Substituting the above definitions of the covariance matrix and the corre- 
lation matrices in equation (2-76) yields the following. 

P n (t) = ' Kn) 'PCt. t-1) P n ^ t_1 ) t-D (I ' KH) T - (I - KH) C (t) D T K T 


- KD C^(t) (I - KH) T + KDW D T K T 


(2-84) 


Substituting (2-82) and (2-83) into (2-84) yields 


P n (t) - (I - KH) (p'(t) + C y (t) 


W" 1 C J(t)j (I - KH) T - (I - KH) C v (t) D T K T 


- KD cj(t) (I - KH) T + KD W D T K T 


(2-85) 


PHIL.CO 
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Using (2-80), the expression C vn (t) W ^ C^(t) may be evaluated 

C vn (t) W " 1 C vn (t) " (I ' C v (t) W ’ 1 C v (t) (I “ ^ 

- (I - KH) C v (t) D T K T - KD <£(t) (I - KH) T 

+ KD W D T K T (2-86) 


Substitution of (2-86) into (2-85) yields 

P n <t) = P’(t) + c vn (t) w‘ l Cvn (t) 


where 


P^(t) - (I - KH) P'(t) (I - KH) T 


(2-87) 


Equation (2-87) is the recursion of (2-81) and shows the feasibility of 
separating the effect of the measurement biases from the effect of random 
errors. The optimal filter gain, K, is to be selected so the trace of 
the covariance matrix, P (t) is a minimum. 


Equation (2-85) can be expanded into the following form. 


p (t) * p'(t) + c w ' 1 c T - k(h[p * ( t) + C 

n v v 1 1 


“- 1 $ 


+ D C 


P'(t) + C 


+ K 


W" 1 C T v ) H 

Mp'(t) + c v w - 1 cl'j 


t + c v d t k t 


T T T 

H + D W D + D C* H 


H C v D T ] 


K 


(2-88) 
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or 


P n (t) » P'(t) + C W* 1 C T • K A • A 1 K 1 + K B K 1 

n ' ' ' ' v v 


(2-89) 


where 


A 

B 


T 

ww 


h(p’(c> c*) 


+ DC 


HP'(t) H + 


|d w + H C v (w' 1 ) T J I 


Dw + H C (w" 1 ) 1 

V 


W. 


Equation (2-89) has the same form as (2-19) with the following equivalent 
variables . 


I ~ P(t) A - H 


K ~ K 


B(t) ~ SS T 


The solution of (2-89) for the optimal gain is given by (2-23). Making 
the appropriate substitutions into (2-23) yields 


K(t) S(t) 


T 

A A (t) 



(2-90) 


or the optimal gain when considering the measurement bias error parameters 
is 


K(t) = A T (t) B _1 (t) 


(2-91) 


where 
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T 

A X (t) 


' ( P ' 


(t) + c w' 1 c, 




C D 
v 




.H P'(t) H + 


[»- + « C v (w-V] 



+ H C (w" 1 ) 
v ' 


f 

— 



-1 


T 

ww = W 


Equation (2-91) is the optimal gain with the measurement bias errors 
considered. 


2.3 PARAMETRIC ERROR ANALYSIS 


The error analysis of a space navigation system generally requires that 
the effects of the error sources on the system performance be examined in 
a parametric manner. A computer simulation of a space mission is in general 
quite complex and the quantity of navigation data to be processed can be 
large. These factors result in a significant amount of computer time being 
required for one data run with a single set of error parameters (variances). 
When there is interest in a large range of variances on each of the error 
parameters this implies a large number of computer runs and the corres- 
pondingly large computer run time. An extremely useful error analysis 
technique would be one by which following a single computer run all the 
parametric data of interest could be obtained simply by means of matrix 
manipulations of the covariance matrix rather than additional runs. This 
section will describe two such parametric analysis techniques that have 
been utilized in this study. The two techniques used apply to: (1) analysis 

of the effects of initial uncertainties in error parameters that have been 
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included in the state and estimated (parameter estimation) and, (2) analysis 
of the effect of errors that have been neglected from the estimation process. 

The parametric analysis techniques being presented were originated by 
Grunckel^ and are elaborated on in Reference 9. 

2.3.1 Initial State Uncertainties 

The change in the state covariance matrix due to the inclusion of an 
observation is given by (2-44). 

P z„ (t) ■ p z (t) - P z (t) ( H z H z p z (t) H z + ») H z p z (t) (2 - 92 > 


The diagonal elements of P (t) define the variances on the vehicle 
state estimate, position and velocity, and any bias errors that are being 
estimated as part of the expanded state vector. 


For the case of a scalar problem, (2-92) becomes 


P 

n 


p .q 

H 2 P o + Q 


(2-93) 


In a conventional error analysis, where one of the elements of P (t) (2-92) 
may be an equation of motion parameter (mass of a body, A.U. etc.) and a 
number of navigation observations are made of the vehicle, the effect of 
changing the initial variances of the bias parameters in the problem would 
be determined by simulating the mission with new parameters. For the 
scalar case, this corresponds to changing the initial variance, P , by 
a factor k and making a computer run. The technique being presented here 
is one in which an additional direct observation of the bias parameter is 
made at the conclusion of the first computer run. 


2-29 


PHILCO 


PHILCO-FQBO CORPORATION 


Space & Re-entry 
Systems Division 



SRS-TR148 


If the initial parameter variance of interest in the scalar case is 


P o = k p 0 (2-94) 

then a direct observation of the bias parameter with a measurement variance, 
Q> defined by equating (2-94) and (2-93) will yield the desired initial 
parameter variance, P 1 . 


P 0 

P P' k P - - ■ ■■ 

n o o P + Q 


(2-95) 


where 


H 1 


Solving for Q yields 


Q 



(2-9 6) 


The fact that the direct observation of the bias parameter with the 
measurement error defined by (2-96) to change an initial variance is not 
required to be made at the start of the run but can be made last is shown 
as follows. 


The order in which the data are processed is shown below, where y (1) is 
the direct parameter measurement. 


PHILCQ 



y(D 

“ H(l) x(l) + q(l) 


y(2) 

* H(2) x(2) + q(2) 


y(t) 

- H(t) x(t) + q(t) 

(2-97) 
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An entirely equivalent set of data would be obtained by processing the 
direct parameter measurement last. 

This gives the following sequence of measurements. 


y(z) = H(2) x(2) + q(2) 
y(t) = H(t) x(t) + q(t) 

y(l) = H(l) cp 1 (t, 1) x(t) + q(l) (2-98) 

The fact that the order of the observations can be interchanged without 

changing the final state uncertainty can be easily seen if the observation 

are all referenced to the same time and the covariance matrix P is 

n 

written as an inverse information matrix (see Appendix A). 

The observation sequence in (2-97) all referenced to time, t, becomes 


- 1 . 


y* (*•) ~ H(l) cp (t, 1) x(t) + q ( 1) = M(l) x(t) + q(l) 
y 2 (t) = H(2) cp -1 (t, 2) x (t) +q(2) = M(2) x(t) + q(2) 


-1 


y t (t) = H(t) cp" (t, t) x(t) + q(t) = M(t) x(t) + q(t) 
Now, using the following inverse form of the covariance matrix 


f 

1 /• ' 


p n (t> - ( ? 0 , 

(t) | 1 + h t q" 1 h'> 


J J 


(2-99) 


( 2 - 100 ) 

and applying it successively to each of the measurements in (2-99) yields 
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P„(t) 


■ 


p Q (t) 


t 

) 

i=1 


M T (i) Q(i) 4 M(i) 


-1 


( 2 - 101 ) 


Clearly, the order of the terms under the summation in (2-101) is not 
important. Since the order of the observations is not important it is 
possible to change an initial variance by means of a direct observation of 
the parameters at the end of a sum. 


When the parameter whose initial variance is to be changed is a bias whose 
value is unknown but constant for the entire mission, the measurement partial 
is 


M(t) H H( 1) c p _1 (t, 1) 



( 2 - 102 ) 


Where it is assumed that the state variable of interest is the i state 
variable . 


Equation (2-92) which shows the computation of the updated covariance matrix 
can be simplified by noting 


P (t) H A (t) = C. 
2 2 1 


v'> H*<t> - <4 


where 


H z (t) s M(t) = 


000 .... 1 00 


2 . 


C. is the ith column of P (t) and cr. is the (i, i) element of P (t) 

X 2 1 z 
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Then using (2-92) 

/ k a 2 \ -1 

P 2n (,:> = • c i K + r^r) c i < 2 - 103 > 

where 

2 j-h 

°xi *" s t * ie uncertainty in the l state variable 

k is the variance scaling desired 

Using (2-103) the initial uncertainties can be varied by orders of magnitude 
with almost no computational penalty. 

A special case of equation (2-103) is with k - 0. In this case the initial 
uncertainty is zero and the state variable is effectively eliminated from 
the estimation process. 


2.3.2 Neglected Bias Errors 


The parametric analysis of the effect of bias errors that have been neglected 
in the estimation process is straight forward. The error sensitivity 
matrices to equation of motion and measurement biases are shown by (2-29) and 
(2-35) respectively. 


C (t) = 

u' 


(I - KH) 


cp(t, t-1) C u (t-1) 


9 u ( fc > t- 


1) 


(2-104) 


C v (t) = (I - KH) cp(t, t-1) C v (t-1) + K D 


(2-105) 


The covariance matrix of the error in estimate due to neglecting these 
errors is the following. 
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P(t) - C u (t) E(uu T ) C*(t) + C v (t) E(w T ) C*(t) (2-106) 


It can be seen in (2-106) that the contribution of each bias error source 
adds algebr iacally into the total error. Each error source can therefore 
be examined separately and its effect analyzed merely by scaling the vari- 
ance of the element up or down. For example, consider the i ^ element in 
the equation of motion errors. 


P.(t) = 


C lu (t) k E^u*) cj u (t) - C iu k a \ <4(t) 


(2-107) 


th 

where C. (t) is the i column of the matrix C (t) and k is the 
ru u 

desired scaling of the variance of the (i, i) element of E(uu T ). 

This type analysis can be done very easily at the end of a run in which 
the C u and C v matrices have been computed. ^ 
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SECTION 3 

ANALYSIS OF BIAS ERROR SOURCES 


This section presents a description of the bias error sources that are 
considered in this study* There are two basic types of bias error sources 
of interest: (1) equation of motion parameters and (2) measurement param- 
eters. These error sources correspond to the system inputs, u and v, that 
are used in the analysis in Section 2. These parameters enter the estima- 
tion process because even though an estimate is available for each parameter 
there is some degree of uncertainty associated with it. The use of the 
term bias error to describe these errors arises from the fact that although 
the expected difference between the parameter estimate and its true value 
is zero, the uncertainty in the estimate causes a constant error. The 
actual difference enters the estimation process as a constant. For example, 
the estimate of the mass of the earth may be described as follows: 

E(u e ) = 8 eI»\ - ujj - £ (3-1) 

£ = estimate of earth's mass 

= variance in the estimate 

2 

The quantity cr^ is a statistical quantity related to the estimation process 
used to obtain £ . If the true value of the earth's mass is k, then the 
planetary mass in a trajectory simulation is in error as follows: 

g = k - £ (3-2) 

k = the true value of mass 
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£ = the estimated value of mass 

e = actual error in estimate 

The error, e, would introduce a constant bias in the simulation using £, 
it would not exhibit a random variation along the trajectory. This example 
describes the nature of the bias errors being considered. They are con- 
stants in the process being modeled that have been estimated but have a non- 
zero uncertainty associated with their estimates. 

The equation of motion error sources that are included in the analysis 
are: (1) astronomical unit conversion, (2) mass of the Earth, (3) mass 

of Mars, and (4) solar radiation pressure. The measurement error sources 
evaluated are a bias in the on-board sextant measurement and a bias in the 
onboard clock. The following sections will describe each of these error 
sources and how they affect the trajectory estimation process. A large 
part of the information presented in this section is taken from Reference 
10, the midterm progress report, entitled, "influence of Uncertainties In 
The Astronomical Unit Conversion and Mars Planetary Mass on Earth-Mars 
Trajectories . " 

3.1 EQUATION OF MOTION BIASES 

The various methods used in estimating the heliocentric and planetary 
constants are described in References 11 through 16. Tables 3-1, 3-2, and 
3-3 summarize the results of several determinations of the ratio of the 
astronomical unit to laboratory units and the planetary masses of Earth 
and Mars . The large discrepancy between the radar measurements of the 
astronomical unit distance and the dynamical method using the asteroid Eros 
is discussed in Reference 14. The discussion indicates that a plausible 
explanation of this discrepancy is the existence of systematic ephemeris 
errors that are not accounted for in the dynamic method. 
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3.1.1 Astronomical Unit Conversion 

The analysis of the astronomical unit (A.U.) presented here and the results 
in Section 5 are an extension and generalization of the work reported in 
reference 17 by S. Henrick, et al. 

The uncertainty in the ratio of the astronomical unit to a laboratory unit 
is an important factor in the accuracy with which an interplanetary mission 
can be performed. Reference 17 demonstrates the importance of using the 
basic "Gaussian" gravitational constant based on the A. U. and the solar 
mass in trajectory computations. This is due to the eight or nine figure 
accuracy' 1 } to which it is presently known. The same constant expressed 
in laboratory units is only accurate to three or four figures. The impor- 
tance of the uncertainty in the ratio to an interplanetary mission is due to 
the fact that with an ephemeris expressed in terms of the A.U., mission 
analysis specifications of injection conditions at Earth are in terms of the 
A.U. The uncertainty in the conversion of the geocentric injection condi- 
tions from a working laboratory unit to the astronomical unit results in the 
Earth escape velocity being in error in units of A.U. /Day. Conversely, the 
uncertainty in the ratio will appear in the initial heliocentric position 
and velocity of the Earth, the gravitational constant, and in the terminal 
position and velocity of Mars, if these quantities are converted from astro- 
nomical units to kilometers. 

The error caused by the ratio uncertainty in the conversion of the tra- 
jectory problem totally into A.U.'s is the same as the error in converting 
the problem to kilometers^ 1 ^ . This equivalence is shown below. The com- 
puter simulation used in the analysis of the uncertainty in the ratio ex- 
presses the problem in kilometers. The geocentric hyperbolic excess 
velocity is assumed to be known precisely and the uncertainty in the ratio 
occurs in the planetary ephemeris. 

The ratio of the A.U. to the equatorial radius is the ”solar parallax” 
expressed in radians (Figure 3-1). Then the desired ratio of the kilometer 
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to the A.U. embodied in the mean Earth distance, R, is related to the solar 
parallax as follows: 


R = 



(3-3) 


where 


a = Earth equatorial radius 

€ 

tt = solar parallax 

R = Earth-Sun mean distance 

The relative uncertainty in the ratio is 


AR = _ Att 

R TT 



(3-4) 


or neglecting the smaller uncertainty in 


a 

e 


AR Att 

R TT 


(3-5) 


The effect of the relative uncertainty, tt' , in the ratio will appear in 
the initial geocentric position and velocity of the vehicle if they are ex- 
pressed in the astronomical unit. The analysis of the error resulting from 
the conversion of the initial state to astronomical units is presented below. 

The Hohmann transfer trajectory geometry is shown in Figure 3-2. The 
relative error in the major axis, 2a, can be found from the vis-viva 
integral which may be written: 
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(3-6) 


where 


r y 


V 


r 


e» 


v 

e 


r 


e 


2 


v 

€ 


2a 


heliocentric position and velocity of vehicle 
heliocentric position and velocity of Earth 

gravitational constant assuming Earth orbit 
is circular 

major axis of transfer 


For the Hohmann transfer, the initial heliocentric vehicle state is the 
following : 


r - r e v = v e * (3-7) 

where is the velocity of the vehicle at about a million kilometers. 

In the process of conversion of the problem from kilometers to A.U.'s 

the position and velocity of the Earth may be assumed to be known accurately 

in astronomical units. 


Ar = Ar g = Av g = 0 (3-8) 

The heliocentric vehicle velocity is in error due to the fact that v 

cn 

although known accurately in laboratory units must be converted to A.U.'s 
using R as a conversion factor. 

Av = Av^ = v m tt' (3-9) 
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Then from equation (3-6), 


A(2a) vAv 

2 2 
(2a) 2 r/ e 


(3-10) 


or using equation (3-9) 


A2 (a) _ 2_a V V qq 

(2a) r 2 

« v e 


(3-11) 


The conversion of the problem from A.U. to kilometers yields the following 
expression*^) for that differs from equation (3-11). 


A2a 

2a 



(3-12) 


This is the uncertainty in Aa expressed in kilometers , leaving the 
position of Mars as an uncertainty. Since the position of Mars is well 
known in astronomical units, the uncertainty should be sought in A(2a/R) 
not A(2a) . Then equation (3-12) becomes 


*(r) 

. (m 

\2a 

ar\ 

. 2a 

V V 

2 rr» 

1 2a, 

R / 

r e 

2 17 

If 

(r J 



V 

e 


(3-13) 


which is in agreement with equation (3-11). 
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Using approximate Hohmann transfer numbers, equations (3-11) and (3-12) 
are evaluated to use as a point of reference for the data that follow. 


2a 

r 

e 


5 

2 



13 

12 


v 

00 

V 

€ 


1_ 

12 


2a 


375- 10 6 km 


2a 

r 

e 


v v 

00 


2 

V 

€ 


.226 


Letting tt' 


.67 -lO -5 (AR 


1000 km) 


From equation (3-11) the uncertainty in the semi-major axis of the transfer 
is the following. 


A(2a) = (.226) (.67- 10 -5 ) (375- 10 6 ) 

A(2a) = 560 km 

A (a) = 280 km 


The uncertainty in the semi-major axis from equation (3-12) which leaves 
the position of Mars as an uncertainty is the following 

A (2a) = (.226-1) ( . 67 • 10 -5 ) (375-10 6 ) 

A(2a) = -1940 km 

A (a) = -970 km 

The Hohmann transfer example case is illustrated in Figure 3-3. The figure 
shows the significance of the "two uncertainties" in the transfer major axis. 
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The computer simulation used to analyze the effect of the A.U. conversion 
uncertainty on the interplanetary trajectory is described in Section 4. 

3. 1.2 Planetary Mass 


The target approach phase of an interplanetary trajectory is a target 
centered hyperbola (Figure 3-4). The characteristics of this trajectory 
are determined by the vehicle velocity state relative to the target at the 
time the sphere of influence 11 ^ ^ is reached and the planetary mass of the 
target body. The vehicle velocity state relative to the target at this 
time is determined by the particular heliocentric transfer trajectory that 
is used. The following is an analysis of the effect of a planetary mass 
uncertainty on an approach trajectory. The analysis follows that presented 
in Reference 20. 

The error in planetary mass is related to an error in the semi-major axis 
of the approach hyperbola through the vis-viva equation. 


where 



(3-14) 


a « semi-major axis 

u = planetary mass 

v^ = hyperbolic excess velocity 

or 

Jfe = Aa (3-13) 

u a 

The angle between the approach and regression asymptotes, 5 , is related 
to the approach trajectory as follows. 
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where 


6=2 cos * ( ) = 2 cos ^ 


, u + r v / 

p OO ' 


(3-16) 


e = eccentricity 


r 

P 


= periapsis radius 


Figure 3-5 presents the scattering angle at Mars as a function of distance 

of closest approach and the hyperbolic excess velocity. Deviations in 6 

due to the uncertainty in the planetary mass may be obtained as follows. 

-* ( 21 ) 

The B vector magnitude is maintained constant i.e.; 


| B | s b = aV |(l-e 2 ) | 


(3-17) 


Then from equations (3-16) and (3-17) 


. .5. ... 2 Ae 

sin (j) A6 — 

e 


(3-18) 


and 


0 = 2a(l-e 2 ) Aa - 2a 2 «Ae 


or 


Ae = * ■- ' £ — ^ Aa = _ 

a e j 


b 2 Aa 


a e 


(3-19) 


Substituting equation (3-18) into (3-19) yields 


A6 = - 


2 b 2 Aa 


(ae) 3 sin j 


. ( b )2 

' a ' l „3 


2 1 Au 

„ ~ ~r — 

e sin j ' 


, £ 2b Au 2 ab Au 

46 * - ~r u 2— r - 

a« a + b 


(3-20) 
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in terms of v^ equation ( 3 - 20 ) becomes 


A6 


2b_ r — 

2 L UT/v + b 

V co 


CO 



(3-21) 


Figure 3-6 presents the deviations in the scattering angle at Mars as a 

function of v ffl and distance of closest approach. The planetary mass un- 
Au 

certainty, — , is 0.0047. These data were obtained by taking the difference 
between a nominal scattering angle and scattering angles obtained when 
using perturbed values of the planetary gravitational constant, 4 , in a 
conic trajectory program. The difference results obtained in this manner 
for the stated planetary mass deviation agree quite well with the linear 
deviation expressed by equation ( 3 - 21 ). 

The importance of these scattering angle deviations on a Earth return trajec- 
tor y i® expressed by the sensitivity of the Earth close approach distance 
to the scattering angle at Mars. This sensitivity for typical Mars-Earth 
trajectories ranges from one hundred thousand to a million kilometers for 
one degree variation in the scattering angle. 

The deviation in Mars close approach distance as a function of planetary 
mass and close approach distance (RCA) is shown in Figure 3-7. The data 
in Figure 3-7 indicate that an uncertainty in the planetary mass of the 
order shown in Table 3-3 causes close approach deviations from ±2 km for a 
high energy trajectory to +15 km for a very low energy trajectory. The 
deviations on the low energy trajectory increase to +35 km for a close 
approach distance of 50,000 km. These data show that the planetary mass 
uncertainty is an important factor for missions requiring terminal accuracies 
on the order of 15 km and less. 

The entry corridor at Mars with a 5 mb atmosphere is approximately 20 kra^ 22 ) 
or +10 km from a nominal trajectory. This indicates that for an atmospheric 
entry mission, a low energy approach trajectory could have significant devia- 
tions due to the uncertainty in the planetary mass. The statistical 
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significance attached to the uncertainties shown in Table 3-3 is also a 
factor in determining the need for approach guidance corrections on the 
higher energy trajectories. If it is assumed that the uncertainties shown 
in the table represent one sigma values; then the deviations in Figures 3-6 
and 3-7 represent the maximum deviations to be expected in 68 % of the cases 
for a selected uncertainty. It would then require the deviation numbers to 
be increased by a factor of 3 to include 99% of the cases. The uncertainties 
in Table 3-3 have been treated as one sigma values in the analysis. 

The deviation in close approach distance at Earth as a function of planetary 
mass deviation and approach energy is shown in Figure 3-8. These data 
indicate that the deviations are less than +1 Km for the planetary mass 
uncertainties shown in Table 3-2. This indicates that for an entry mission 
at Earth the planetary mass uncertainty is not a very significant factor. 

The data shown in Reference 23 indicates that the entry corridor at Earth 
is 21 Km for a vehicle with a L/D ratio of one and a speed of 18 km/sec. 

3.1.3 Solar Radiation Pressure 

In the trajectory analysis of a spacecraft, a knowledge of the environmental 
forces acting on the vehicle is required. An important component of the 
spacecraft environmental force is that due to radiation. While radiation 
is generally considered as a mode of energy transfer, it is well known that 
momentum is also transferred by radiation. The forces resulting from this 
momentum transfer can have an important influence on the trajectory of a 
spacecraft . 

Several papers have appeared dealing with the effects of solar radiation 
pressure in regard to space vehicles^ 2 ^ ,2 “^ . The theoretical basis of 
radiation pressure and methods for exact determination of the force acting 
on a body in a stream of radiation are discussed in Reference 26. Figure 
3-9 shows the solar radiation pressure variation as a function of distance 
from the sun. 
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The primary interest in this study centers on the effects of uncertainties 
in the solar radiation force acting on the vehicle. The analysis presented 
in Section 5 treats the solar radiation force uncertainty as a bias error* 
The analysis examines the ensemble behavior of the vehicle state through 
the state convariance matrix. This analysis requires the use of the 
control sensitivity matrix, ^ LL ( t 2 > t 1 )> as described in Section 2. This 
matrix is obtained by integration of the variational equations. These 
equations are described in Section 3.1.4. The solar pressure enters the 
equations of motion as follows. 


The equation of motion of the spacecraft may be written as follows: 


where 


X = F (X , U, t) 

r 

R 
V 

Using the partitioned form of equation (3-22) 


(3-22) 


R = F'(R, V, U, t) 


(3-23) 


The function F* may be written as follows: 

R. 


R N / A i 

F ' ■ - 1 * ? ■ 1=1 u * {7 + 


/U. K. \ K 

(?*?)* 


(3-24) 
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where 

r is the magnitude of R, the vector from the central body 
to the vehicle. 

is the gravitational mass of the central body. 

M' i is the gravitational mass of the i th body. 

T ± is the magnitude of R^ , the vector from the central body 
to the i th body. 

6^ is the magnitude of the radius vector from the i^” 
body to the vehicle. 

Pj are peturbation accelerations due to drag, oblateness, solar 
pressure, etc. 


The forms of the solar radiation pressure in the p. term of (3-24) is 
assumed to be an inverse square repulsive acceleration relative to the sun. 


P 


- k -f 
sp 3 


(3-25) 


where 

R g is the radius vector from the sun to the vehicle. 

r is the magnitude of R . 
s s 

k gp is a constant including the reflectance of the vehicle, 
the vehicle's mean sun-directed surface area and the 
vehicle's mass. 


An uncertainty in the constant, k gp , of (3-25) produces an uncertainty in 
the solar radiation pressure acceleration. 
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3.1.4 Variational Equation 

The statistical analyses for both the solar radiation pressure and planetary 
mass uncertainties presented in Section 5 make use of the control sensi- 
tivity matrix ^ u ( t 2 > • This section presents a description of the varia- 
tional equations that are integrated to obtain both and 

In order to find the sensitivity of the state at time t 2 to the state devia- 
tions and variations in the equation of motion parameters at t^, the varia- 
tional equations are integrated. 


The acceleration of the vehicle was given in (3-23) as 

R = F'(R, V, U, t) 


Any small variation in the V, R, or U result in a first order variation in 
R given by 

®‘IF® + w :sv + lr 4u <3 ‘ 26 > 

where means the gradient of F' with respect to the vector A. 


The partial derivative of R with respect to any arbitrary parameter, ot, is 
therefore : 


BR _ *F' 3R 3V ?5F ' jHJ 


(3-27) 


by the chain rule of differentiation. If the order of differentiation in 
equation (3-27) is interchanged, the following second order linear differen- 
tial equation results. 


d 2 , 3R x _ $F 1 . SFj_ d_.dR. . 3F' Su 

, 2^' ^R '■^R ; dt''^ + 5U ?o> 

at 


(3-28) 
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Now if Of represents the initial state, then -gj = 0 and the solution of 
(3-28) subject to the following initial vehicle state conditions: 



t-o 


(3-29) 


provides the state transition matrix, $(t,o) * % ■ / y . The sensitivity 

9x(t) °x\°) 

of the state to a control, ^(t,o) = — — L y which is meaningful only if 
the control vector is constant over the time interval under consideration, 
is obtained by solving equation (3-28) with ot = u, subject to the following 
initial conditions. 



t-o 





(n - dimension of control vector) 


(3-30) 


3.2 MEASUREMENT BIASES 

The measurement biases considered in the analysis of Section 5 are: (1) 

onboard clock bias, and (2) an angle bias in the sextant measurements. 
Section 2 shows that the manner in which these bias errors enter the 
estimation process is through the gradient of the measurement with respect 
to the biases, D(t). The following is a derivation of the measurement 
gradients for the sextant type planet-star observations that are used in 
the analysis. The measurement sensitivities to the following quantities 
are derived. 
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R Position state-radius vector from the central body 
to the vehicle. 

V Velocity state-inertial velocity relative to the central body 
t Time bias-error in onboard clock 
6 Angle bias-error in sextant instrument 

The sextant star-planet angle measurement is made in the plane defined by 
the vector from the vehicle to the star and the vector from the vehicle 
to the body being observed. Figure 3-10 illustrates this measurement. 

The selection of suitable stars to use in making these measurements is 
discussed in Reference 27. 


From Figure 3-10 the sextant measurement is the following: 


O' = cos S • P) = cos l ( a T #) 


(3-31) 


where 


P = 


The total measurement gradient that is desired is the following: 


. do> dcr dcr ifrr 
z " ( dR ’ ; *5 ; ' 


(3-32) 


Taking the gradient of (3-31) with respect to R yields 


• ^ o T ^ ^ 

sin a l* = S Hr + ? HU 


(3-33) 


Since 


3R 
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Equation (3-33) becomes 


to 


** |p j sin o 


Ci - p £ T ] = 


J — A A 

1 r _L^_s_ x f 
p L |p * II 


The gradient of (3-31) with respect to V yields 


(to 

¥ 


= o 


The gradient of (3-31) with respect to the bias angle 6 is 


■53 " 1 


The gradient of the measurement with respect to the time bias is 


to to dP 
St ~ ^ ^ 


or 


to S 

Tto = <- 7 


£-s> (i7T Cl ' * ' v ^) 


Using (3-35) equation (3-39) can be written as 


Tr “ |[< P p * V > 


The total measurement gradient with respect to the state and the 
is shown below. 


H * ' {' 


pI 


j xj 

| A A 

UPxS 


x P 


; 0 ; 1 ; 


*0 


PxS 
|PxS | 


[p - v] 
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(3-35) 


(3-36) 


(3-37) 


(3-38) 


(3-39) 


(3-40) 
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SECTION 4 

DIGITAL COMPUTER PROGRAMS 


The analysis of the bias errors described in Section 3 has been performed 
using three digital computer program simulations. The "Patched Conic A.U. 
Program" and "Planetary Mass Program" are simulations of the trajectory 
in two dimensional conic form. The analysis performed using these programs 
is in the form of deviations or differences from a nominal trajectory due 
to deviations in the physical constants. The third program used is the 
"Mark II Error Propagation Program."^®) This program is designed for use 
in making a statistical error analysis of a navigation system that utilizes 
a Kalman filter in the data processing. 

The following sections will describe the capabilities of these three programs. 
The Mark II description is limited to those features that are used in the 
study . 

4.1 PATCHED CONIC A.U. PROGRAM 

The Patched Conic A.U. Program simulates the effect of a change in the A.U. 
Conversion to a laboratory unit on the planetary ephemeris and shows the 
resulting effect on an interplanetary trajectory. 

The planet ephemeris model used in the program has the following character- 
istics. The planets Earth and Mars are on coplanar circular orbits about the 
Sun at distances of 1 A.U. and 1.53 A.U. respectively. The uncertainty in 
the A.U. conversion* is included in the model in the following manner. For 
the Earth on a two body Keplerian orbit, the A.U., mass of the Sun, mass of 
the Earth, and the period of the Earth about the Sun are related by the 
following expression. 


*The nominal conversion factor used is 149599000 km/A.U. 
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u , 2 - 


G (M + M ) 
N a e_ 

(AU) 3 


(4-1) 


where 


U) 

M 

s 

M c 

AU 


= angular frequency of the Earth 
= mass of the Sun 
= mass of the Earth 
= astronomincal unit 
~ universal gravitational constant 


It is assumed that the Earth's angular frequency is known perfectly and 
that the earth s mass can be neglected with respect to the Sun's mass. 
Under these assumptions, the partial derivative of equation (4-1) becomes 


z' GM 

^ AU J = 3( AU“ ) 

UJ = CONST u 


(4-2) 


The relationship shown in equation (4-2) indicates that a change in the 
"length" of the A.U. must be accompanied by a change in the mass of the 
Sun in order to maintain U) constant. In the ephemeris model used, a change 
in the A.U. is accompanied by changes in the radial distances of the planets 
and the mass of the Sun. These changes maintain the angular frequencies 
of the planets constant. 


The interplanetary trajectory generation is accomplished as follows. The 
launch and target planets are positioned with an initial angular separation 
that will satisfy the planetary geometry required for a specified heliocentric 
transfer angle and flight time. 

The trajectory program obtains an Earth-Mars (Mars-Earth) heliocentric conic 
trajectory with a specified flight time and transfer angle. The heliocentric 
conic is then patched to a Mars (Earth) centered conic trajectory at the 
sphere of influence. The initial heliocentric velocity magnitude is then 
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varied in a differential correction loop to obtain a specified close 
approach distance at Mars (Earth) . This process establishes a nominal 
trajectory for the flight time and transfer angle. The initial heliocentric 
velocity vector is then separated into two parts as shown below. 


v = 



+ v 


00 


where ^ 

v = initial heliocentric vehicle velocity 
— * 

Vg = heliocentric velocity of Earth 

v oo = geocentric hyperbolic excess velocity 

(velocity relative to Earth at a million km) 


(4-3) 


The geocentric hyperbolic excess velocity, v o , represents the Earth depar- 
ture condition measured in kilometers/sec. This is the initial heliocentric 
velocity that a mission analysis would show is required for a nominal 
ephemeris. This velocity is assumed to be known precisely and is not changed. 
The A.U. conversion factor is then perturbed causing changes in the positions 
and velocities of Earth and Mars. The gravitational constant is also changed 
in accordance with equation (4-2) . The result of these changes is that the 
initial vehicle state relative to the Sun deviates from the nominal conditions 
The vehicle position is changed with the change in the Earth's position. 

The vehicle velocity relative to the Sun is changed through the change in 
the Earth's velocity in equation (4-3). The perturbed heliocentric tra- 
jectory is patched to Mars and the approach trajectory differences from the 
nominal computed. The process is shown pictorially in Figure 4-1. 


This program can call the guidance program described in Section 4.4. The 
velocity required to correct the approach trajectory differences are then 
computed as a function of time. 
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4.2 PLANETARY MASS PROGRAM 

The Planetary Mass Program is designed to analyze the effect of uncertainties 
in planetary mass on interplanetary approach trajectories. The character- 
istics of a planetary approach trajectory are determine by the vehicle 
velocity state relative to the planet at the "sphere of influence" and 
the planetary mass. 

The trajectory model used in the program is a conic section. During the 
approach phase of a mission, this is a good approximation to the three 
dimensional trajectory. 

The B vector and radius of closest approach (RCA) are used to describe the 

vehicle passage of the planet. The B vector and the associated unit vectors 

ft, ft, f (Figure 4-2) are described in Reference 21. The S vector is in 

the direction of the approach asymptote and the R, f vectors are in the 

. £ “♦ 

plane normal to the S vector and containing the B vector. 

The magnitude of the vehicle velocity state at the sphere of influence, v^, 
is a program input used to simulate approach trajectories of different 
energies. The flight path angle of the approach velocity vector is used 
as a control by the program to obtain a trajectory with a specified distance 
of closest approach. 

After the generation of the desired nominal approach trajectory, the program 
perturbes the planetary mass and computes the differences between the per- 
turbed approach trajectory and the nominal. The program can call the 
guidance program which will compute the velocity required to correct the 
differences . 
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4.3 MARK II ERROR PROPAGATION PROGRAM 

The Mark II Error Propagation Program is used in the statistical analyses 
of the effects of the bias errors on the navigation system performance. 

The capabilities of the program that are used in the analysis and the 
program modifications added during the study are described below. 

The program has the capability of simulating an onboard navigation system 
with the following types of measurements: 

a. Range (Radar) 

b. Range Rate 

c. Theodolite (Right Ascension and Declination) 

d. Sextant (Star-Planet Angle) 

e. Range (Subtended Angle) 

The sextant measurement is the only one used in the study. The star-planet 
angles used in the analysis are measured in two specific directions relative 
to the trajectory plane. The angles being measured are assumed to lie in 
the trajectory plane or normal to it. Figure 4-3 illustrates the orienta- 
tion of the two angles and shows the stars required for such measurements . 

The standard deviation, <?, of the instrument random error is input to the 
program. The instrument noise is the quantity, q, in the expression for a 
measurement, y, shown below. 

y = Hx + q (4-4) 


where 

y measurement 
H measurement gradient 
x Deviation state 
q random noise 
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The program is capable of analyzing the effects of the following bias errors. 

a. Sextant Bias (Modification added during study) 

b. Onboard Clock Bias 

c. Earth Planetary Mass Uncertainty 

d. Mars Planetary Mass Uncertainty 

e. Solar Radiation Pressure Uncertainty 

The standard deviation on the estimate of each of these quantities is input 
to the program. 

The Mark II program, with modifications incorporated during the study, is 
capable of analyzing the effects of the bias errors for each of the three 
techniques described in Section 2.2; (1) Neglecting, (2) Including, and 
(3) Considering. The parametric analysis techniques described in Section 
2.3 have also been added to the program capability. 

The error analysis quantities computed by the program and used in Section 5, 
are defined below with a summary of the equations used in the Kalman trajec- 
tory estimation and end point prediction processes. A basic assumption in 
the theory is that linearity is satisfied in the neighborhood of a nominal 
trajectory. 

Between the onboard observations, the deviation state estimate and the 
error covariance matrix are propagated in time along the nominal trajectory 
as follows. 


x(t 2 ) = tPCtjjtp x( tl ) 


(4-5) 


P(t 2 ) = cp(t 2 ;t 1 )P(t 1 )cp T (t 2 ;t 1 ) 


(4-6) 
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At the time of an observation, the measurement information is included 
in the state estimate and a new covariance matrix obtained in the following 
manner . 

*n = *o + K <y-y> (4-7) 


where 


P „ * p „ - ““'o < 4 - 8 > 

K = p 0 ^ T (^ p 0 hT + Q) (Kalman Filter Gain) 

H = measurement gradient 
y - actual measurement 

A 

y ~ estimate of measurement 

Q = convariance matrix of random measurement noise 


The equations shown above are for the case in which all the bias errors are 
neglected. The sequence of program operations described by equations (4-5) 
through (4-8) is the same when using the different types of treating bias 
errors but the form of each equation is different (see Section 2). 


The Mark II program computes the statistical quantities at any requested 
time point along the nominal trajectory. The program also computes the 
statistical characteristics of the estimate of the trajectory end point 
constraints. The estimated deviation state and covariance matrix of the 
error in estimate are propagated to the end point. The propagation to the 
end time, T, appropriately includes the effects of bias errors in the 
equations of motion if they are being analyzed in the simulation. 


x(T) = cp( T ,t) x(t) (4-9) 

P(T) = 9(T,t) P(t) cp T (T,t) (4-10) 

where 

T = End point time 
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The estimated end point deviation, x(T), and the covariance matrix of the 
error in estimate, P(T) , are then transformed into an estimate of the 
constraint deviations and the associated error in estimate respectively. 

The constraints used^to describe the end deviations in the program are 

B T and B R. The constraint data are obtained by means of a point 

transformation as follows. 


6 B • T 

I - G(T) x(T) 

6 B • R, 


(4-11) 


/ B-T pa 1 CT 2 \ 

Vr ) ' G(T) P(I) cT(t > <4 - 12) 


where 


G(T) 


b B-T \ 
b x(T) 

S B-R 

a x(T)/ 


point transformation 


Equations (4-5) through (4-12) describe the processes by which the statis 
tical navigation data shown in Section 5 are obtained. 


4.4 GUIDANCE PROGRAM 


The guidance program is used by the A.U. Program and the Planetary Mars 
Program to compute the velocity required to correct specified trajectory 
deviations. The program is capable of using two guidance laws; (1) fixed 
time of arrival (FTA) , and (2) variable time of arrival (VTA). The three 
end point constraints used with each guidance law are shown below: 


PHILCO BUS 



1 X(T)\ 


1 b-t' 

FTA = 

Y(T) 

VTA - 

B.R 


\ Z(T)/ 

T 
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i# 


where 

T - nominal arrival time 

“* A A 

B*T,B*R = orthogonal components of the B vector 
= hyperbolic excess velocity 
X,Y,Z = nominal vehicle position state at time, X 

The guidance velocity correction required at time, t, is computed in the 
following manner. The vehicle deviation state is propagated along a nominal 
trajectory to the end point (equation 4-13) and transformed into appropriate 
constraint deviations (equation 4-14). 


x(T) = <p(T,t) x(t) 


D(T) = C(T) x(T) = C(T) cp(T,t) x(t) 


(4-13) 

(4-14) 


where 

— • 

x = deviation state vector 

D = constraint deviation vector 

C(T)= point transformation from the state to either 
FTA or VTA constraints 

The sensitivity of the end constraints to a velocity correction at time, t, 
is obtained from the partioned transition matrix. 


I 

I 

» 

V 

I 

I 


A(T, t) = (A ; A ) = C(T) cp(T,t) 


(4-15) 


3x6 


3x3 3x3 3x6 6x6 


where 


sensitivity of end constraints to a position change at 
time, t. 

sensitivity of end constraints to velocity change at 
t ime , t . 
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The velocity correction required to null the constraint deviation vector, 
D(T) , in equation (4-14) is the following. 

— ♦ 

A 2 * g <t) + D(T) = 0 
or 

-4 

^g (t) * A 2 L = -( A 2 1a ! x > *<*> (4-16) 

— * 

where x the guidance velocity correction. 

O 

The deviation state, x(t), used in the guidance program is input by the 
parent program. It is obtained by taking the difference between a nominal 
approach trajectory and a perturbed trajectory. The trajectory is perturbed 
due to an uncertainty in the planetary mass or an uncertainty in the astro- 
nomincal unit conversion. 
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SECTION 5 

MIDCOURSE NAVIGATION STUDY 

The results presented in this section show the effects of bias errors on the 
performance of an onboard navigation system. These results are an evaluation 
of the navigation system requirements generated in the original study^ with 
the constraint that only random errors were considered in the analysis. The 

analysis is performed for the midcourse phases of both the outbound Earth 
to Mars trajectories and the return Mars to Earth trajectories. 

The navigation analysis is limited to the use on onboard sextant measuring 
device with a random error of 10 arc seconds. The initial study showed 
that an instrument of approximately this accuracy was capable of performing 
the Earth Mars round trip navigation such that ±3.5 km entry corridors were 
achieved at both Mars and Earth. The onboard measurement schedules used 
for the outbound and return trajectories are shown in Tables 5-1 and 5-2 
respectively. The analysis of the schedule selection is in Reference 1. 

There are several trajectories utilized in the following analyses. The 
statistical analyses with the Mark II program uses the same round trip 
nominal trajectory used in the original study. This nominal trajectory, 
a 235 day outboard and 296 day return with a 40 day stay at Mars, is shown 
and described in Appendix B. The trajectory was selected from the study 
results described in Reference 29. The Patched Conic A.U. Program and 
Planetary Mars program are used to analyze equation of motion error effects 
on a large number of trajectories. The trajectories analyzed include many 
transfer and approach trajectories of practical interest. 

The initial covariance matrices of the error in estimate used in the statis- 
tical analyses on the outbound and return trajectories are shown in Tables 
5-3 and 5-4 respectively. 
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5.1 EARTH-MARS TRANSFER 

The results of two forms of analysis are presented below for the outbound 
trajectories. The first type of analysis evaluates the trajectory devia- 
tions or differences from a nominal trajectory that are introduced by 
deviations in the equation of motion parameters. These analyses are per- 
formed with the Patched Conic A.U. Program and the Planetary Mass Program. 

The second type of analysis is a statistical evaluation of the effect of 
bias errors on the navigation system performance. These results are obtained 
with the Mark II Error Propagation Program. 

5.1.1 Equation of Motion Errors 

The equation of motion bias errors analyzed are uncertainties in the 
following: (1) A.U. conversion, (2) Earth and Mars planetary masses, and 

(3) solar radiation pressure. 

5.1. 1.1 A.U. Conversion . The data and results presented in this section 
were generated using an uncertainty in the conversion of the A.U. to kilo- 
meters of ±1000 km. This is slightly larger than the uncertainty shown in 
Table 3-1 for the 1963 adopted value. 

The approach phase of a number of Earth-Mars trajectories is analyzed to 
determine the navigation and guidance requirements due to the uncertainty 
in the A.U. conversion. Five heliocentric transfer angles are used with 
flight times for each from 100 to 500 days. Three trajectories of interest 
listed in Table 5-5, are included in the analysis. They are: (1) Ho hma nn 

transfer 180 degrees, 260 days, (2) Mariner IV trajectory 160 degrees, 228 
days, and (3) High energy outbound leg of round-trip trajectory 270 
degrees, 235 days. 

The trajectory data and corresponding target approach deviations are shown 
in Figures 5-1A through 5-5B. Part A of each figure shows the flight time, 
launch velocity, and target approach velocity as a function of the direction 
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of the hyperbolic approach asymptote, S. Part B shows the deviation in 
close approach distance for a 1000 km change in the A.U. conversion to 
kilometers. The 160, 180, and 200 degree transfers each show two deviation 
minimums. The 225 and 270 degree transfers each has a single minimum. 

The minimum deviations are near zero for the 180 degree transfer and increase 
with transfer angles away from 180. The minimum deviations for the 225 and 
270 degree transfers are 125 km and 550 km respectively. These data show the 
possibility of selecting trajectories that minimize the effect of the 
uncertainty in the A.U. conversion on the close approach distance. The 
Mariner IV trajectory is one that is near a minimum. The 228 day 160 degree 
transfer has a deviation of 175 km uncertainty in the conversion. 

The position deviation state at the sphere of influence (patch point) for all 
the trajectories is approximately 1000 km. The close approach deviation 
minimums are the result of these errors at the patch point being in direc- 
tions that result in cancellation or partial concellation of the deviation 
in the periaries distance. 

The trajectories marked with an asterisk on the 160, 180, and 270 degree 
transfers are analyzed to determine the approach guidance velocity required 
to correct the deviations. The results of this analysis are shown in 
Figures 5-6, 5-7, and 5-8. The solid lines indicate the Av required for a 
fixed time of arrival (FTA) and the dotted lines the requirements for a 
variable time of arrival (VTA). The curves show that the requirements for 
FTA are nearly the same for all the trajectories shown. A correction at 
the sphere of influence requires about 10 meters/second and grows to 
approximately 200 meters/second as periaries is approached. The tsv require- 
ments for the VTA guidance law show a wide variation depending on the specific 
trajectory selected. 

The trajectories that have small close approach deviations, curve 1 in 
Figure 5-6 and curves 1 and 3 in Figure 5-7, have velocity requirements 
that range from less than 1 meter/second at the sphere of influence to 8, 4 
and 2 meters/second respectively at periaries. The remaining VTA curves in 
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Figures 5-6 and 5-7 show larger Av requirements that range from 2 meters/ 
second to 40 meters/second for the trajectories with larger deviations. 

These requirements are considerably smaller than those required with a FTA 
guidance law. Figure 5-8 shows the Av requirements on a 270 degree transfer. 
Figure 5-5B shows the minimum deviation for the 270 degree transfer is 
550 km, which is much larger than other transfer minimums and the velocity 
requirements are correspondingly higher. The VTA velocity requirements 
are only slightly smaller than those required for a FTA guidance law. 

The time at which a reasonable guidance correction can be made is determined 
by the navigation system. The error in estina te of the end constraints 
must be below a predetermined level before the guidance maneuver can be 
executed. The selection of an entry mission at Mars with a ±10 km entry 
corridor defines tolerable limits on the end constraint deviations. If it 
is required that the confidence in hitting the entry corridor is to be 
99 percent (3 sigma), then the one sigma error in estimate of the close 
approach distance must be reduced to ±3.3 km. The guidance correction can 
then be made with a 99 percent confidence (neglecting execution errors) of 
hitting the ±10 km corridor. 


Figure 5-9 shows the error in estimate of the end constraints for three 
approach trajectories of different energies. The navigation measurements 
are made with a 10 arc second sextant. Measurements are taken every 15 min 
minutes. The initial error in estimate of state is assumed to be 1000 km 
in each of the inplane position coordinates and 0.2 meters/second in the 
velocity coordinates. These errors correspond to the actual deviations 
that occur at the time of patch to the target due to a 1000 km uncertainty 
in the A.U. conversion. Due to the onboard observations, the error in 
estimate of the constraints is quickly reduced to less than 100 km. It then 
remains relatively constant until the last few hours of the approach. The 
error in estimate is sufficiently small for an entry mission (3.3 km) 
approximately 3 or 4 hours prior to periaries on each trajectory. 
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Figures 5-6 through 5-8 indicate that this time corresponds to corrections 
of 50 to 70 meters/second for a FTA guidance law. The VTA guidance require- 
ments at this time are less than 10 meters/second except for the 270 degree 
transfer where they are about 30 meters/second. A FTA guidance policy 
allowing for two approach corrections could reduce the total Av required 
considerably from the 50-70 meters/second required for a single correction. 

A factor that has been neglected in the guidance analysis is the execution 
errors. A correction of 70 meters/second with proportional errors of 1 
percent would produce a 0.7 meter/second execution error. Three or four 
hours before periaries the close approach sensitivity to a velocity change 
is such that 0.7 meter/second error will cause deviations that are the same 
order of magnitude as the entry corridor. This factor also favors the 
guidance policy of two smaller approach corrections for an accurate planet 
passage. 

5. 1.1.2 Mars 1 Planetary Mass . The results presented in this section show 
the navigation and guidance requirements for controlling the approach 
trajectory under the influence of an uncertainty in the planetary mass. 

The results assume that the midcourse guidance system has controlled the 
vehicle to the sphere of influence perfectly. The only equation of motion 
uncertainty considered is the planetary mass. 

The time history of the growth in the predicted deviations in close approach 
distance and B magnitude (equations 4-9 and 4-11) based on the state devia- 
tion is shown in Figure 5-10. These data were obtained using a planetary 
mass uncertainty of 130 km /sec and a close approach distance of 5000 km. 
The curves all display the characteristic of having very small deviations 
until 4 to 8 hours before periaries. The deviations then grow rapidly to 
values from 2 to 15 kilometers. The approach guidance Av required to 
correct these deviations is shown in Figure 5-11 as a function of time along 
the trajectory. The requirements are shown for both FTA and VTA guidance 
laws. The Av required on these trajectories for each guidance law is 
between 1 and 10 meters/second during the last few hours. The VTA velocity 
requirements are smaller in all cases. 
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The time at which a guidance correction can be made is determined by the 
capability of the navigation system to estimate the end point deviations 
to a satisfactory accuracy. The capability of an onboard navigation system 
using a 10 arc second sextant to estimate the end constraints is shown in 
Figure 5-12. The results are shown for three nominal trajectories with 
different energies. The parameters being estimated include the vehicle 
state and the planetary mass. The initial vehicle state uncertainty is 
assumed to be zero and the uncertainty in the planetary mass is 130 km /sec. 

The tolerable error in estimate for an entry mission is shown on Figure 5-11 
as +3.3 km. The times on these trajectories that this level is reached 
are 2 days 19 hours for the trajectory with = 2.0 km/sec and 1 day 12 hours 
for the - 4.0 km/sec trajectory. Using these correction times in Figure 5-11 

shows the guidance velocity requirements are approximately 1 meter/sec for a VTA 
guidance law and 3 meters/sec for FTA guidance law. The significance of these 
approach corrections in terms of the total mission is discussed in Section 6. 

5. 1.1.3 Statistical Analysis . The statistical analysis of the navigation 
system is performed with the Mark II Error Propagation Program. The measure- 
ment instrument is a 10 arc second accuracy sextant. The results presented 
in this section are obtained while using the initial state covariance matrix 
shown in Table 5-3 and the measurement schedule shown in Table 5-1. The 
sequence of measurements used within the schedule is the following. Three 
inplane measurements are made followed by a single out of plane measure- 
ment. This sequence is continually repeated. The nominal trajectory used 
is described in Appendix B. 

The results presented in Figures 5-13A and B show the effects of neglecting 
the equation of motion bias errors on the error in estimate of the end point 
constraints B*T and B-T respectively. The T vector is in the trajectory 
plane and the R vector is normal to the trajectory plane. The theory for 
the error analysis being used is presented in Section 2.2.1. 

The number 4 curve in each figure shows the expected error in estimate of 
the end constraint under the assumption of a perfect physical model. These 
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uncertainties shown are due to injection errors and random errors in the 

measurements. These results correspond to those obtained in the original 

( 1 ) 

study. The three remaining curves in each figure show the additional 

error in estimate due to the neglected uncertainties in the planetary 

masses of Earth and Mars and the solar radiation pressure. The planetary 

3 2 3 2 

mass uncertainties used are 15 km /sec for the Earth and 150 Km /sec for 
Mars. These are approximately the uncertainties shown in Tables 3-2 and 
3-3 for the 1961 adopted values. The assumed form of the solar radiation 
pressure acceleration from equation (3-25) is 


P 


k 

sp 



(5-1) 


The solar radiation constant, k g p, * s a function of the vehicle mass, 
vehicle area projected normal to the vehicle-sun line, and the surface 
reflectance. For a hypothetical vehicle with a reflectance of one, a 
constant normal area of 100 square feet, and a weight of 200 pounds, the 


constant, k 


sp 


7 3 2 

is approximately 10 km /sec . 


The results shown in 


Figures 5-13A and B are for an uncertainty in the solar radiation pressure 
c o o 

constant of 10 km /sec or 10%. 


The characteristics of the three curves in Figures 5-13A and B showing 
the error due to neglecting the uncertainties in the planetary masses and 
solar radiation pressure are due to the following effects. The initial 
uncertainties show the error each would introduce if no navigation measure- 
ments were taken. At different points in time, each of the errors begins 
decreasing due to the measurements that are being taken. The decrease in 
error begins to occur at a time when the equation of motion error source 
has had some influence on the trajectory and the navigation system is esti- 
mating the perturbed trajectory. For example, the planetary mass of Earth 
has an immediate strong effect on the trajectory. A large portion of the 
effect of this error is quickly removed by the navigation system. In 
contrast, the mass of Mars has essentially no effect on the trajectory until 
the sphere of influence is reached (234 days). Therefore, the errors in 
the end constraint estimates due to neglecting this error remains constant 
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(number 2 curves) until the last day. At this time the mass of Mars does 
influence the trajectory and the effects of its uncertainties can be 
removed . 

The total error in estimate of the end constraints due to neglecting the 
three equation of motion parameters is shown in Figures 5-14A and B by 
the number 1 curves. These curves are obtained by taking the root sum 
square of the errors due to the individual error sources. Curve number 2 
in Figures 5-14A and B is the error in estimate when the bias errors are 
neglected. Curve Number 3 in these Figures is the total error in estimate 
of the constraint. This curve is the root sum square of curves 1 and 2. 

The difference between curves 2 and 3 shows the degree to which the results 
of the original study concerning the navigation system performance 
were optimistic due to neglecting the bias errors. 

The B*T constraint estimate is the more critical of the two because it 
indicates the accuracy to which the entry altitude is known. As indicated 
earlier, for an entry mission at Mars the altitude corridor is *3.5 km. 

“ ♦ A 

Figures 5-13A and 5-14A show that the error in B’T due to neglecting the 
uncertainty in Mars planetary mass is 4 to 5 km. The effect of an error 
of this magnitude must be evaluated in terms of mission requirements. 

While this error makes an entry mission at Mars marginal, it would have 
only a minor effect on a flyby or orbiter mission. 

During the approach phase of the mission, Figures 5-13A and B indicate that 
the effects of neglecting uncertainties in the planetary mass of Earth 
and the solar radiation pressure are negligible compared to the effect 
of Mars' mass uncertainty. This result indicates that these two error 
sources can be neglected without affecting the end point estimation per- 
formance. 
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Bias Error Included . The equation of motion bias errors instead of being 
neglected can be included in an expanded state vector and estimated in 
addition to the vehicle state. The theory describing this estimation 
process is presented in Section 2.2.2. The results obtained with this 
type of processing are shown in Figure 5-15A and B and 5-16A and B. 

Figures 5-15A and B show the effects of an uncertainty in the solar radia- 
tion pressure when it is being estimated as part of the state. Parametric 

curves are shown for uncertainties in the solar radiation constant ranging 

6 3 2 ^ 

from 0 to 10 km /sec . The curve for a zero uncertainty in the constant 

is equivalent to neglecting the bias error. The differences between this 
zero uncertainty curve and the curves for the solved for bias error show 
the extent ot which the end point estimate is optimistic when neglecting 
the radiation pressure uncertainty. These differences approach zero at 
approximately 230 days. At this time, the error in estimate of the end 
constraint for the zero uncertainty case is no longer optimistic. Essen- 
tially the same conclusion can be reached from the data shown in Figures 
5-13A and B. 

The error in estimate in Figures 5-13A and B due to neglecting the solar 
pressure uncertainty completely is of significant magnitude for a 
slightly longer time (234 days) and then becomes negligible. This is 
due to the fact the uncertainty is not being estimated as is the case for 
the results shown in Figures 5-15A and B. 

The effects of uncertainties in the planetary masses of Earth and Mars 
are shown in Figures 5-16A and B. These results were obtained while 
including the mass uncertainties in the state vector and solving for them. 

The conclusions that can be drawn from these results are in general agree- 
ment with those shown in Figure 5-13A and B for the neglected mass uncertain- 
ties. The effect of uncertainties in the Earth's mass is negligible after 
approximately 200 days. At this time, the curves for the zero uncertainty 
and the maximum uncertainty form a single curve (Figures 5-16A and B) . The 
second half of the curves in Figures 5-16A and B show the effect of Mars 
planetary mass uncertainty. The error in estimate of B’T in Figure 5-16A 
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shows that uncertainties from 0 and 450 km^/sec^ in the planetary mass 

cause the error in estimate to increase from 3 km to approximately 20 km. 

— * ^ 

An error as large as ±20 km in the B-T estimate would be unsatisfactory for 
an entry mission at Mars. 

5.1.2 Measurement Bias Errors 

The analysis of the effects of the measurement biases is performed with the 
Mark II Error Propagation Program. The bias errors are analyzed by two 
techniques. One technique is that of including them as part of the state 
and estimating them. This analysis technique is described in Section 2.2.2. 
The second technique used is that of considering the effect of a bias error. 
This process is described in Section 2.2.3. The Mark II program is not 
capable of analyzing the error due to completely neglecting measurement 
biases as was done for the equation of motion errors. 

The results shown in Figures 5-17A and B are for the two methods of treating 
the bias error in the sextant measurement. Curve one shows the error in 
estimate of the end constraint for an instrument with a 10 arc second 
random error and no bias error. Curve two shows the error in estimate for 
the addition of a 10 arc second bias that is included in the state and 
estimate. The third curve is for the case in which the 10 arc second bias 
is considered in tne estimation process but is not estimated (variance is 
held constant) as part of the state. Curves one and two in Figures 5-17A 
and B show that the increase in the error due to an instrument bias error 
that is solved for is negligibly small (less than 1 km at the end time). 

Figure 18 shows the manner in which the standard deviation of the measure- 
ment bias is reduced as it is estimated along with the state. The figure 
shows the reduction for initial uncertainties of 10 and 20 arc seconds. 

In both cases, the initial uncertainty is quickly reduced to two arc seconds 
and finally to one arc second. 
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Curve three shows that there is a considerable increase in the error in 
estimate when the bias error is considered but not estimated. The end 
point estimate error in B-f is increased from 7 km to 14 km. This result 
indicates that if the error source were neglected completely, the error 
due to neglecting it would be larger than the 7 km difference shown above. 

Figure 5-19 shows the effect of a 60 second onboard clock bias on the error 
in estimate. The clock bias is included in the state and estimate. These 
results indicate that the clock bias can be estimated rapidly and its effect 
removed from the estimation process. The effect of completely neglecting 
this error source is not presented. The high correlation between the clock 
and the estimate, as shown by the ability to solve for the bias, indicates 
it would be a significant error if present and neglected. 

5.2 MARS -EARTH TRANSFER 

The results obtained on the effects of bias errors on the Mars-Earth 
trajectories are presented in this section. The data presentation is 
restricted because the results and conclusions that can be reached are in 
general the same as obtained for the outbound portion of the mission. 

The analysis of the measurement bias errors is omitted. The analysis in 
Section 5.1.2 indicates that these are error sources that can not be neglected 
without introducing significant errors. This conclusion applies to the 
return trajectories equally well. 

5.2.1 Equation of Motion Errors 

The equation of motion bias errors analyzed are uncertainties in the 
following: (1) A.U. Conversion, (2) Earth and Mars planetary masses, 

and (3) solar radiation pressure. 
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5.2. 1.1 A.U. Conversion . The data presented in this section were computed 
using a conversion uncertainty in the A.U. to kilometer of ±1000 km. This 
is slightly larger than the uncertainty shown in Table 3-1 for the 1963 
adopted value. These data were generated with the Patched Conic A.U. 

Program. 

Five heliocentric transfer angles are used with flight times for each from 
100 to 500 days. One trajectory of interest included in the analysis is 
the Mars-Earth nominal trajectory (296 days 270 degrees) used in the original 
study. The trajectory data and corresponding target approach deviations 
are shown in Figure 5-20A through 5-24B. Part A of each figure shows the 
flight time, launch velocity, and target approach velocity as a function 
of the direction of the hyperbolic approach asymptote, S. Part B shows 
the deviation in close approach distance for a 1000 km change in the A.U. 
Conversion to kilometers. The 160 and 180 degree transfers each show two 
deviation minimums. The remainder of the curves have a single minimum. 

These are essentially the same type of characteristics that where exhibited 
on the outbound trajectories shown in Figures 5- 1A through 5-5B. 

The 270 degree transfers show a very narrow deviation minimum around the 
point of a 470 day flight time. The 296 day nominal trajectory from the 
first study has a close approach deviation of over 2000 km for the 1000 km 
conversion deviation. 

The conversion of these trajectory deviations into velocity requirements 
was not performed. The velocity requirements, as in the cases for the 
outbound trajectory, would vary from near zero at the minimum points to 
extremely large values away from the minimums. 

5.2. 1.2 Statistical Analysis . The statistical analysis of the navigation 
system is performed with the Mark II Error Propagation Program. The 
measurement instrument is a 10 arc second sextant. The results presented 

in this section are obtained while using the initial state covariance matrix 
shown in Table 5-4 and the measurement schedule shown in Table 5-2. The 
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sequence of measurement within the schedule is three inplane measurements 
followed by a single out-of-plane measurement. The nominal trajectory used 
is described in Appendix B. 

The results presented in Figures 5-25A and B show the effect of neglecting 
the equation of motion bias errors on the estimate of the end constraints 
BT and B-R respectively. The theory of the error analysis being used is 
presented in Section 2.2.1. 

The number 4 curve in each figure shows the expected error in estimate of 

the constraint under the assumption of a perfect physical model. These 

uncertainties are due to the injection errors and random errors in the 

measurements. These results correspond to those obtained in the original 

study. ^ ^ The three remaining curves in each figure show the additional 

error in estimate due to the neglected uncertainties in planetary masses 

of Earth and Mars and the solar radiation pressure. The planetary mass 

uncertainties used are 15 km 3 /sec 2 for the Earth and 150 km 3 /sec 2 for 

Mars. These are approximately the uncertainties shown in Tables 3-2 and 

3-3 for the 1961 adopted values. The uncertainty in the solar radiation 

constant, k , is taken as 10 6 km 3 /sec 2 . 
sp 

The curves in Figures 5-25A and B have the same character as those shown 
in 5-13A and B for the outbound trajectory with the roles of Earth and 
Mars planetary masses interchanged. In the return case, the uncertainty 
in the mass of Earth is so small that it has negligible effect on the 
terminal accuracy. The mass of Mars and the solar radiation pressure 
exhibit an effect until late in the flight. During the final day, their 
effects become quite small. 

The total error in estimate of the end constraints due to neglecting the 
three equation of motion parameters is shown in Figures 5-26A and B by 
the number 1 curves. These curves are obtained by taking the root sum 
square of the errors due to the individual error sources. Curve number 2 
in Figures 5-26A and B is the error in estimate when the bias errors are 
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neglected. Curve number 3 in these figures is the total error in estimate 
of the constraint. This is the root sum square of curves 1 and 2. The 
difference between curves 2 and 3 shows the degree to which the results of 
the original study were optimistic concerning the navigation system per- 
formance when neglecting the bias errors. The curve differences indicate 
that neglecting the equation of motion biases would not adversely effect 
an entry mission requiring a ±3.5 km terminal accuracy. 
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SECTION 6 

SUMMARY AMD CONCLUSIONS 

The theoretical analysis in Section 2 presents some new developments in 
the analysis of the effects of bias errors on an orbit estimation process. 
These developments concern three areas of interest. They are the 
following: (1) the effect of neglecting bias errors in the modeling of 

the physical process, (2) the separable properties of the effects due to 
random errors and those due to bias errors, and (3) techniques for 
efficient parametric analysis by means of matrix manipulations. 

The data results obtained on the effects of equation of motion bias error 
sources on the navigation system performance indicated the following. The 
uncertainty in Mars planetary mass produces deviations in close approach 
of 10 to 20 km for practical Mars approach trajectories . The error in 
the state estimate due to neglecting this uncertainty is on the order of 
4 to 5 km at the end point. This is a significant error when considering 
an entry mission with a 3.5 km corridor requirement. The effect of Mars 
mass uncertainty on the return trajectory is negligible following the 
navigation measurements that are used. 

The uncertainty in the planetary mass of Earth causes approach deviations 
of less than 1 km on the practical approach trajectories. On the Earth- 
Mars outbound trajectory, the effects of this mass uncertainty is removed 
by the navigation measurements . 

The effects of a solar radiation pressure uncertainty of 10 percent can be 
removed on both the outbound and return trajectories by means of the navi- 
gation measurements . 

The effect of an uncertainty in the A.U. conversion to laboratory units is 
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a strong function of the particular heliocentric trajectory being used. 
The analysis of 5 heliocentric transfer angles for various flight times 
for both the outbound and return trajectories shows one or two minimums 
in the close approach deviations for each transfer angle. The deviation 
minimums vary from near zero to 550 km for a 1000 km uncertainty in the 
A.U. The minimum deviations are near zero for a 180 degree transfer and 
increase for larger and smaller transfer angles. The deviations can be 
estimated by the navigation system to a satisfactory accuracy for an 
entry mission. 

The guidance analysis is restricted to the determinations of the velocity 
required to correct the deviations caused by the Mars planetary mass 
uncertainty and the A.U. conversion uncertainty. The tSJ required is a 
function of the time at which the correction is applied. The CU required 
to control the deviations due to Mars mass uncertainty for a FTA 
guidance law vary from 5 to 10 meters/second and from 1 to 5 meters/ 
second with a VTA guidance law. The velocity requirements for an 
uncertainty in the A.U. are quite trajectory dependent. The guidance 
velocity corrections for a FTA guidance law are from 50 to 70 meters/ 
second when using only one correction. The corrections for a VTA law 
vary considerably. The trajectories with small deviations (less than 100 
km) require corrections from 1 to 10 me ter s /second . The trajectories 
with the larger deviations require corrections of 10 to 30 meters /second. 

The guidance requirements for an Earth-Mars mission obtained in the 
original study while neglecting the two uncertainties that have been 
described above are shown in Table 6-1. The results in Table 6-1 for a 
VTA guidance law include the effects of errors in an onboard navigation 
system and guidance system execution errors. The approach trajectory 
deviations due to a planetary mass uncertainty cannot be estimated until 
the last few hours of the approach trajectory. It would therefore be 
necessary to control these deviations with the final correction. The 1 
meter/second final correction shown in Table 6-1 would increase to a 
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maximum of approximately 5 meters/second with a mass uncertainty of 
3 2 

150 km /sec . The trajectory deviations due to the uncertainty in the 
A.U. conversion can be estimated with an error of 30 to 40 km one day 
prior to periaries with a 10 arc second instrument. This allows the 
possibility of making a correction at this time that will correct the 
deviations to an accuracy consistent with the estimate. The deviations 
remaining after the correction could then be removed with the final 
maneuver. On the trajectories with large deviations due to the A.U. 
conversion, this would increase the third correction of Table 6-1 by 
approximately 10 meters /second. The final correction would be increased 
by 2 to 3 meters/second. 

The discussion of guidance Av requirements above is summarized in Table 6-2. 
These results were obtained by algebraically adding the velocity 
requirements caused by the two uncertainties in the equation of motion to 
those due to injection errors, navigation errors, and guidance system 
execution errors. This very pessimistic analysis of adding these 
independent effects algebraically increases the total velocity 
requirements from 23 meters/second to 41 me ters /second . 

The sextant angle measurement bias and onboard clock bias both have a 
significant influence on the navigation system accuracy. If these errors 
are neglected, the results imply a degradation in the end point estimate 
accuracy of greater than 7 km. This is not satisfactory for an entry 
mission. If the bias error sources are included as part of the state 
being estimated their effect can be eliminated by the "calibration" of 
the sextant instrument and clock. 
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SECTION 7 
RECOMMENDATIONS 

Additional study areas that would extend the scope of the present study 
and are considered important for defining the navigation and guidance 
requirements of an interplanetary mission are the following: 

a * Guidance Analysis . The guidance analysis performed with 

trajectory deviations in this study should be extended to a 
statistical analysis that includes the effects of equation of 
motion bias errors . 

b. Measurement Biases . The effect of neglecting the measurement 
bias errors should be analyzed. The theory is described in 
section 2.0 but the digital simulation program is not available 
at present. 

c - A. U. Conversion . The theory required for a statistical analysis 
of the effect of the uncertainty in the A. U. conversion should 
be developed. The theory should account for the changing of 
the reference body centers. 

d. Filtering Techniques . The study has assumed the use of a Kalman 
filter in the data processing. The use of other filtering tech- 
niques should be evaluated and, in particular, consideration 
should be given to their onboard implementation. The performance 
of some non-optimum data processing filters used with an 
imperfect physical model should be investigated. 

e * Onboard Computers . The design of the onboard computer should be 
studied to determine means of trading off speed for reliability 
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(500 - 600 day missions). Also techniques for simplifying 
calculations for estimating and predicting the state should be 
investigated. This investigation should include the effects of 
truncation errors in the computer. 

f. Beacons . The importance of having beacons on Mars should be 
evaluated for the approach phase of the mission, terminal 
maneuvering phase and orbital phase . 

g- Powered Flight . The guidance system inertial equipment 
requirements should be determined for the retro maneuver 
(powered and/or atmospheric) and the powered flight out of Mars 
orbit. A retro analysis is also required at perigee on the 
return. 

h. Mars Orbit . The navigation requirements in orbit should be 

determined in terms of specific mission objectives. The influence 
of the oblateness of Mars on these requirements should be evaluated. 
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APPENDIX A 

EQUIVALENT FORM FOR THE COVARIANCE MATRIX 

The state covariance matrix following an observation that is weighted with 
the optimum gain, K(t) is the following. 

P R = P - PH T (HPH T + Q)' 1 HP (A- 1) 


This can be written as 


P n = P - PH T Q _1 (I + HPH T Q -1 ) -1 HP (A- 2) 

Using the matrix identity 

(I + M)" 1 = I - M + M 2 - M 3 + (A-3) 


Equation (A-2) can be expanded 

P n = P - ph t q _1 •[ I - hph t q _1 + [hph t q -1 ] 2 

- [hPH T Q - 1 ] 3 + j HP (A-4) 

f T - 1 T-1T-1 T - 1 1 1 

P n " P \ 1 ' H Q HP + H Q HPH Q HP - (H Q HP) + - (A-5) 

Again using the identity of (A-3) 

P n * P[l + H T Q _1 HP] -1 (A-6) 

or 

P n = C P' 1 + hV^]" 1 (A-7) 

Equation (A-7) shows the equivalence between the Kalman filter and the weighted 
least sequences. 
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APPENDIX B 

INTERPLANETARY TRAJECTORIES 

This appendix presents a description of the two nominal interplane- 
tary trajectories which have been used in this study. The trajectories 
are: (1) Nominal high-energy Earth-Mars, (2) Direct-return Mars-Earth. 
The data which are shown are in Earth equator and equinox of 1950 coor- 
dinates. Ecliptic projections of the trajectories are shown in Figures 
B-l and B-2. 


B.l NOMINAL HIGH-ENERGY EARTH-MARS TRAJECTORY 


This trajectory (Figure B-l) has a launch date of 10 February 1975, 
1 hours, 32 minutes, 28.629 seconds, with a park orbit length of 
1736.518 seconds. 


Date 

Fractional Date 
x (km) 

y 

z 

x (km/sec) 

y 

z 


Earth-Centered 

Conic 

(Injection) 

10 February 1975 
2 h l m 25. 147 S 
-0.51940523+04 
-0.33714096+04 
-0.21758862+04 
0.97623321+01 
-0.11540529+02 
-0.5422576+01 


Sun-Centered 

Conic 

(Patch) 

10 February 1975 
23 h 43. m '24.812 S 
-0.11441701+09 
0.84261958+08 
0. 36545250+08 
-0.95000336+01 
-0.27582250+02 
-0.1185094+02 


Mars -Centered 
Conic 
(Patch) 

10 October 1975 
21 h 35 m 48.596 S 
-0.44850600+06 
0.27226599+06 
0.20978100+06 
0.43146753+01 
-0.26728354+01 
-0.20487671+01 

Mars of 3860 km 
near the ecliptic 


The trajectory has a radius of closest approach at 
(500 km altitude). It passes Mars on the Sun light side 
plane. The flight time is 235 days. 
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B.2 DIRECT-RETURN MARS -EARTH TRAJECTORY 


This trajectory (Figure B-2) is the return trajectory which leaves 
Mars 40 days after the arrival of the nominal high-energy trajectory. 


Date 

Fractional Date 
x (km) 

y 

z 

x (km/sec) 

z 

z 


Mars-Centered 

Conic 

(Injection) 

12 November 1975 

10 h 47 m 26.241 S 

-0.34075300+04 

-0.18867350+04 

-0.66730232+03 

0.32159301+01 

-0.48195511+01 

-0.32001703+01 


Sun-Centered 

Conic 

(Patch) 

13 November 1975 


0.91568038+08 

0.18887866+09 

0.84164875+08 

-0.17725844+02 

0.80037536+01 

0.3498172+01 


Earth-Centered 

Conic 

(Patch) 

3 September 1976 


-0.88066700+06 

0.27005599+06 

0.84061749+05 

0.99308927+01 

-0.30949005+01 

-0.10445975+01 


19 h 41 m 18.707 S 5 h 4 m 57.839 S 


The trajectory has a perigee radius of 6442 km (76 km altitude) . 
Perigee is at a latitude of -59.98 degrees and longtitude of 223.06 
degrees. The flight time is 296 days 18 hours. 
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FIGURE 2-1 PHYSICAL PROCESS 
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FIGURE 2-3 PHYSICAL PROCESS 
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FIGURE 3-6 MARS' SCATTERING ANGLE DEVIATION 
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FIGURE 3-7 DEVIATIONS IN CLOSEST APPROACH DISTANCE 






FIGURE 3-8 DEVIATIONS IN CLOSEST APPROACH DISTANCE 
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FIGURE 3-10 


SEXTANT MEASUREMENT GEOMETRY 
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FIGURE 4-1 EPHEMERIS GEOMETRY CHANGE WITH A.U. CHANGE 










FIGURE 5-IA TRAJECTORY CHARACTERISTICS - 160 TRANSFER, EARTH-MARS 

















180 TRANSFER, EARTH -MARS 







FIGURE 5-3A TRAJECTORY CHARACTERISTICS - 200 TRANSFER, EARTH-MARS 







FIGURE 5-3B DEVIATION IN CLOSE APPROACH - 200 TRANSFER, EARTH-MARS 











FIGURE 5-4 B DEVIATION IN CLOSE APPROACH - 225 TRANSFER, EARTH-MARS 






FIGURE 5-5A TRAJECTORY CHARACTERISTICS - 270 TRANSFER , EARTH-MARS 
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FIGURE 5-5B DEVIATIONS IN CLOSE APPROACH - 270 TRANSFER, EARTH-MARS 





























RMS ERROR IN ESTIMATE OF END CONSTRAINT (KM) 
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FIGURE 5-9 ERROR IN ESTIMATE O! 
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figure s-ii av required for approach guidance 
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FIGURE 5-I6B EFFECT OF PLANETARY MASSES ON B • R ESTIMATE 
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FIGURE 5-20A TRAJECTORY CHARACTERISTICS - 160 TRANSFER, MARS-EARTH 








FIGURE 5-20B DEVIATIONS IN CLOSE APPROACH - 160 TRANSFER, MARS - EARTH 
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FIGURE 5-2 IA TRAJECTORY CHARACTERISTICS - 180 TRANSFER, MARS-EARTH 








FIGURE 5-21 B DEVIATIONS IN CLOSE APPROACH - 180 TRANSFER, MARS-EARTH 




FIGURE 5-22A TRAJECTORY CHARACTERISTICS - 200 TRANSFER, MARS" EARTH 























FIGURE 5-24B DEVIATIONS IN CLOSE APPROACH -270 TRANSFER, MARS-EARTH 
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mars-earth trajectory^ 

RANDOM SEXTANT 


ERROR IN ARC SECONDS 


measurement 


SCHEDULE TABLE 5-2 


BIAS errors 


CD 150 KM 3 / 


Q) 15 KM 3 / 


© I - 10 6 km 3 /sec 


® biases neglected 


solar pressure (D 


(T) MARS MASS 


( 2 ) EARTH’S mass curve less 
THAN I KM 


IRMHI 






TIME FROM INJECTION (DAYS) 


figure 5-2SA 

error in estimate of b-t due to neglecting bias 

errors 
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FIGURE 5-26B COMPOSITE EFFECT OF BIAS ERRORS ON B-R ESTIMATE 
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Figure B-l Ecliptic Projection, Earth-Mars Trajectory 
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Figure B-2 Ecliptic Projection, Mars-Earth Trajectory 
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TABLE 3-1 


MARS MASS (SUN'S MASS - 1) 


Mars Mass (m~ l ) 

Method 

Author 

Reference 

3648000 

Vesta 

Leveau (1890) 

(3) 

3093500 


Newcomb (1895) 

(3) 

3601280 

Vesta 

Leveau (1907) 

(3) 

3085000 + 5000 

Weighted Mean 

de Sitter (1938) 

(3) 

3110000 + 7700 

Eros 

Rabe (1949) 

(3) 

3079000 + 5700 

Deimos 

Urey (1952) 

(3) 

3090000 + 10000 


Adopted (1961) 

(3) 

3088000 + 3000 

Weighted Mean 

Clemence (1961) 

(8) 

3090000 + 3000 


Adopted (1963) 

(8) 


TABLE 3-2 

EARTH + MOCA) MASS (Sun's Mass - 1) 


Mass (m *) 

Method 

Author 

Reference 

329390 


Newcomb (1895) 

(3) 

327900*200 

Weighted Mean 

de Sitter (1938) 

(3) 

328390*103 

Eros 

Witt (1933) 

(3) 

328452*43 

Eros 

Rabe (1949) 

(3) 

328446*43 

Eros Revised 

E. Rabe (1954) 

(3) 

328440*40 


de Vaucouleurs 
(1961) 

(3) 

328450*50 


Adopted (1961) 

(3) 

328905.2*5 


Adopted (1963) 

(8) 
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TABLE 3-3 
ASTRONOMICAL UNIT 


A. U. (KM) 

Method 

Author 

Reference 

149662400 + 25600 

Eros (Geometric) 

(1941) Jones 

(4) 

149530300 + 10200 

Eros (Dynamical) 

(1950) Rabe 

(4) 

149598640 + 250 

Venus (Radar) 

(1962) Muhleman et.al. 

(5) 

149597850 ± 400 

Venus (Radar) 

(1962) Pettingill et.al. 

(5) 

149598100 + 400 

Venus (Radar) 

(1962) Muhleman (Revi- 
sion of Pettingill's 
Value) 

(5) 

149601000 + 5000 

Venus (Radar) 

(1961) Thompson, et .al. (GB) 

(5) 

149599500 + 800 

Venus (Radar) 

(1961) Kotelnikov (USSR) 

(5) 

149599244 + 278 

Mariner II Tracking 

(1963) Anderson, et.al . 

(7) 

149599000 + 700 

Recommended Value 

(1963) 

(8) 
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TABLE 5-1 
SCHEDULE 


TIME 

BODY 

OBSERVED 

INTERVAL 
OF OBS. 

0-3 h 

Earth 

15 Minutes 

3 h -l D 

Earth 

1 Hour 

1 D -90 D 

Sun 

5 Days 

90 D -180 D 

Mar 8 

5 Days 

180 D -230 D 

Mars 

1 Day 

230 D -233 D 

Mars 

6 Hours 

233 D -234 D 

Mars 

1 Hour 

234 D -234 D 18 h 

Mars 

30 Minutes 

234 D -18 h -234 D 20 h 

None 

— 

234 D 20 h -234 D 22 h 

Mars 

15 Minutes 

234 D 22 h -235 D 

None 

— 


TOTAL NUMBE 
OF OBS. 



TOTAL 
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TABLE 5-2 


SCHEDULE 


SRS-TR148 


Body 


OBS Interval 


Number of 
OBS 



l d -40 d 
40 d -242 d 
242 d -292 d 
292 d -295 d 
295 d -296 d 
296 d -296 d 12 h 
296 d 12 h -296 d 14 h 
296 d 14 h -296 d 16 h 
296 d 16 h -296 d 18 h 


Earth 


Earth 


Earth 


Earth 


Earth 


Earth 


Earth 


15 Minutes 


1 Hour 


2 Days 

5 Days 
1 Day 

6 Hours 
1 Hour 


30 Minutes 


15 Minutes 




Total Number of OBS 
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TABLE 5-4 

INJECTION COVARIANCE MATRIX OP STATE DEVIATIONS PROM NOMINAL 


° n = 17.3 Km, a v = 17.3 Km, a * = 17.3 Km, 

C 1 = 17.3 m/sec , a- =17.3 m/sec, <*• = 17.3 m/sec 
n v w 


n v w 


a 


v 


w 


n 

v 

w 

n 

v 

w 


Symmetric 

Matrix 


0 

1 


0 0 
0 0 
1 0 
1 


0 

0 

0 

0 

1 


0 

0 

0 

0 

0 

1 


TABLE 5-5 

TYPICAL EARTH-MARS TRAJECTORIES 


Transfer 

Trajectory 


Flight Time 

(Day 8 ) 


Heliocentric 
Angle (Deg) 


Mars Approach 
Velocity (KM/SEC) 


Hohmann 
Mariner IV 
Outbound of 
Round Trip 


260 

228 

235 


180 

160 

270 


2.6 

3.1 

6.6 
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TABLE 6-1 

GUIDANCE PERFORMANCE VTA GUIDANCE LAW 


Correction 

End Constraint 
Deviations 
(KM) 

Av Req'd 
M/Sec 

Time 


-* * 
B-T 

-* A 

B-R 



1 

10300 

2390 

10.56 

1 D 

2 

269 

153 

8.18 

220° 

3 

12.8 

8.6 

3.53 

234 D 

4 

6.54 

2.01 

.92 

234°20 h 


TABLE 6-2 

VELOCITY REQUIREMENTS WITH EQUATION 
OF MOTION UNCERTAINTIES 


Correction 

Av Req'd 
M/Sec 

i 

10.56 

2 

8.18 

3 

3.53 

10.00 (AU) 

4 

0.92 

5.00 (U) 

3.00 (AU) 

TOTAL 

41.19 
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